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ABSTRACT 

We present results from the Spitzer /IRS spectral mapping observations of 15 local luminous infrared 
galaxies (LIRGs). In this paper we investigate the spatial variations of the mid-IR emission which in- 
cludes: fine structure lines, molecular hydrogen lines, polycyclic aromatic features (PAHs), continuum 
emission and the 9.7 /im silicate feature. We also compare the nuclear and integrated spectra. We find 
that the star formation takes place in extended regions (several kpc) as probed by the PAH emission 
as well as the [Ne II] 12.81 /im and [Ne III] 15.56 /im emissions. The behavior of the integrated PAH 
emission and 9.7 jam silicate feature is similar to that of local starburst galaxies. We also find that the 
minima of the [Ne III] 15.56 /im/[Ne II] 12.81 /im ratio tends to be located at the nuclei and its value is 
lower than that of H ll regions in our LIRGs and nearby galaxies. It is likely that increased densities 
in the nuclei of LIRGs are responsible for the smaller nuclear [Ne III] 15.56 /im/[Ne II] 12.81 /im ratios. 
This includes the possibility that some of the most massive stars in the nuclei are still embedded 
in ultracompact Hll regions. In a large fraction of our sample the 11.3 /im PAH emission appears 
more extended than the dust 5.5 jam continuum emission. We find a dependency of the 11.3 /im PAH/ 
7.7 /im PAH and [Ne II] 12.81 /im/11.3 jam PAH ratios with the age of the stellar populations. Smaller 
and larger ratios respectively indicate recent star formation. The estimated warm (300 K < T < 1000 
K) molecular hydrogen masses are of the order of lO^M©, which are similar to those found in ULIRGs, 
local starbursts and Seyfert galaxies. Finally we find that the [Ne II] 12.81 /im velocity fields for most 
of the LIRGs in our sample are compatible with a rotating disk at ~kpc scales, and they are in a 
good agreement with Ha velocity fields. 

Subject headings: galaxies: nuclei — galaxies: structure — galaxies: active — galaxies: starburst — 
infrared: galaxies 



1. INTRODUCTION 

Infrared (IR) br ight galaxies were first identified al- 
most 40 years ago (jRieke Lowlll972f ). The /i?^ 5 satel- 
lite detected a larg e number of these galax ies in the 
local Universe (see iSanders fc Mirabel 1 19961 for a re- 
view). Galaxies were classified according to their IR 
(8—1000 /im) luminosities into Luminous Infrared Galax- 
ies (LIRGs, IO^^Lq < LiR < lO^^Lo) and Ultra Lumi- 
nous Infrared Galaxies (ULIRGs, Ljr > lO^^L©). Since 
then, an increasing number of studies showed the impor- 
tant role of (U)LIRGs at cosmological distances. LIRGs 
and ULIRGs are not common in the local Universe, ac- 
counting for just 5% and < 1%, respectively, of the to- 
tal IR emission of galaxies. However, at z~l, LIRGs 

'pereira^damir. iem.csic.es 

* This work is based on observations made with the Spitzer 
Space Telescope, which is operated by the Jet Propulsion Labo- 
ratory, California Institute of Technology under NASA contract 
1407 

^ Instituto de Estructura de la Materia, CSIC, Serrano 121, 
E-28006, Madrid, Spain 

^ Associate Astronomer, Steward Observatory, University of 
Arizona, AZ 85721, USA 

^ Steward Observatory, University of Arizona, 933 N. Cherry 
Avenue, Tucson, AZ 85721, USA 

^ Actual address: University of Crete, Department of Physics, 
GR-71003, Heraklion, Greece 

^ Ritter Astrophysical Research Center, University of Toledo, 
2801 West Bancroft Street, Toledo, OH 43606, USA 

^ Departamento de Astrofisica, Facultad de CC. Fisicas, Uni- 
versidad Complutense de Madrid, E-28040 Madrid, Spain 



dominate the IR background and the co-moving star 
formation rate densit y, wh ile ULIRGs are dominant at 
z^2 (ILe Fl oc^h et al.l [2QQ5I : iPerez-Gonzalez et al.l [20051 : 
iCaputi et a l. 2007). 

Thanks to the high sensitivity of the Infrared Spec- 
trograph (IRS) (.Houck et al. 2004) on-board Spitzer^ we 
can, for the first time, study systematically the mid- 
IR properties of (U)LIRGs, both locally and at high- 
z. The majority of IRS local studies hav e, so far, fo- 
cused o n the pro perties of starburst galaxies (Smith e t al. l 
2007k iBrandTet a l. 2006) and ULIRGs (Armus e tHI 
120071 : iFarrah et aDl2007[ ). The latter class of galaxies 
includes a large variety of physical and excitation condi- 
tions. For instance, most of these local {z < 0.3) ULIRGs 
are powered predominantly by compact starburst events, 
although a substantiaLAGN contribution is present in 
about half of theni pmanishi et al.l 120071 : IFarrah et all 
120071 : iNardini et alll2008f ). In the future the larger sam- 
ple (^200) of local LIRGs of the Spitzer legacy program 
GOALS will enable a statistical study of this class of 
galaxies (jArmus "eFall 12009). 

Although ULIRGs at z~2 extend to even higher lumi- 
nosities than local ones, their mid-IR spectra are similar 
to those of local starbursts, rather tha n to local ULIRGs 
(jRigbv et al.ll2008l : IFarrah et al.ll2008D . A si milar behav- 
ior is found for high-z submillimeter gal axies (|Pope et al.l 
120071 : iMenendez-Delmestre et al.l l2009f ). This behavior 
may indicate that star formation in IR-bright galaxies 
at z~2 is taking place over larger scales (a few kpc) 
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than in local ULIRGs, where most of the IR emission 
arises fr om very compact regions (sub-kilopa rsec scales, 
see e.g., iSoifer et al.ll2QQQl : iFarrah et al"]l2QQ7[ ). 

As a result, local LIRGs may provide important in- 
sights to the behavior of high-redshift ULIRGs. Their 
star formation is often d istributed over a substantia l frac- 
tion of the galaxy (e.g.. lAlonso- Herrero et al.ll2QQ6[ ). and 
in any case is usually not embedded at the high optical 
depths characteristic of the nuclear star forming regions 
in local ULIRGs. To make good use of this analogy re- 
quires observations that map the local LIRGs to get an 
understanding of their global properties as well as to pro- 
vide a better comparison with the observations of high-z 
galaxies, where the entire galaxy is encompassed in the 
IRS slit. The spectral mapping mode of Spitzer /IRS is 
well suited for this purpose, since it can provide both the 
nuclear and integrated spectra of local galaxies, enabling 
study of the spatial distribution of the spectral features 
that make up the integrated spectrum. 

This is the first in a series of papers studying the mid- 
IR properties of local LIRGs. We present Spitzer /IRS 
spectral mapping observations of a representative sample 
of fifteen local LIRGs. The first results and the goals 
of this program were presented by I Alonso- Herrero et al.l 
(|2Q1Q| ). In the present paper we focus on the analysis 
of the spatially resolved measurements. The paper is 
organized as follows. In Section [2] we present the sample 
and the observational details. Section [3] describes the 
analysis of the data. We examine the silicate feature in 
Section HI Sections O [6l [71 discuss the spatially resolved 
properties of the atomic fine structure lines, PAHs, and 
molecular hydrogen emission, respectively. We explore 
the velocity fields in Section [H Finally, in Section [9] we 
summarize the main results. Throughout this paper we 
assume a flat cosmology with Hq = 70 km s~^Mpc~^, 
= 0.3 and = 0.7. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. The sample 

The galaxies studied here were s elected from the sam- 
ple of 30 local LIRGs described in lAlonso-Herrero et al.l 
(|2006l ). which was draw n from the IRAS Rev ised Bright 
Galax y Sample (RBGS, [S anders e t al.|[2QQ3l ). The sam- 
ple of [Alonso- Herrero et a l. (2006) is a volume limited 
sample, and the velocity range was selected so that 
the Paa emission line could be observed with the NIC- 
MOS F190N filter on-board the Hubble Space Telescope 
(HST). This restricted the sample to local LIRGs (40 
Mpc ^ d < 75 Mpc). From this sample we selected 15 
galaxies (Table [T]) with extended emission in the NIC- 
MOS Paa images to be mapped with the IRS spectro- 
graph. This selection criterion does not introduce any 
bias in the luminosity, distance and nuclear activity dis- 
tributions of the mapped galaxies as c ompa red to the 
parent sample of lAlonso- Herrero et al.[ ([20 06'). The re- 
maining galaxies in the sample have been observed using 
the IRS staring mode as part of different programs. The 
latter will be analyzed with the mapping data in a forth- 
coming paper (Pereira-Santaella et al., in preparation). 
The interacting galaxy Arp 2 99 was studied in detail by 
[Alonso-Herrero et al.[ ()2009b[ ) and will be included in the 
appropriate comparisons and figures in this paper. 

Most of the nuclei of the LIRGs in our sample are 
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classified as H ll region-like from optical spectroscopy, 3 
galaxies have a Seyfert classification and 2 are classified 
as compo site (interm ediate between H ll and LINER, see 
lAlonso- Herrero et al.ll2009a) . The infrared luminosity of 
the sample ranges from 10^^"^ to 10^^"^ L©. The main 
properties of the galaxies are summarized in Table [H 

2.2. Spectral Mapping 

We obtained Spitzer IRS spectral mapping observa- 
tions through two guaranteed time observation (GTO) 
programs P30577 and P40479 (PI: G. H. Rieke). We used 
all four IRS modules: high resolution (R~600) modules 
with Short-High (SH; 9.9 - 19.6 /im) and Long-High (LH; 
18.7 - 37.2 /im), and low resolution (R^60 - 120) with 
Short-Low (SLl; 7.4 - 14.5 /im, SL2; 5.2 - 7.7 /im) and 
Long-Low (LLl; 19.5 - 38.0 /im, LL2; 14.0 - 21.3 /im). 
Maps were made by moving the telescope pointing posi- 
tion perpendicular to the long axis of the slit with a step 
size half of the slit width. The SL, LL and LH slits are 
long enough to cover the full region of interest with just 
one slit length, but due to the smaller size of SH slit, 
parallel steps were also needed. The parallel step size 
was half of the slit height. This strategy produced maps 
with a redundancy per pixel of 2 for SL, LL and LH, and 
4 for SH. Table [T] details the integration time, pixel size 
and step size for each IRS module. The maps cover at 
least the central 20" x20" to 30" x30" of the galaxies. 

Dedicated background observations are not needed for 
the low resolution modules (SL and LL) because the sep- 
aration on the sky between the two orders of each module 
is larger than the galaxy size. That is, when order one 
is on the galaxy, order two is observing the background 
and vice versa. The SH and LH slits are too small to 
allow background subtraction. We obtained dedicated 
background observations for LH by observing a region 
2 arcminutes away from the galaxy. SH sky observa- 
tions are available only for the galaxies observed under 
the program P40479. For the galaxies without dedicated 
SH background observations, we estimated the sky back- 
ground emission by combining the background spectra 
from the low resolution modules. Generally the sky con- 
tribution is not very important for the galaxies in our 
sample, except for NGC 3110 and Zw 049.057. The sky 
background is bright at the position of these galaxies and 
its subtraction is essential to match the fluxes from the 
different slits. Moreover, subtracting background obser- 
vations is always useful because it removes many of the 
rogue pixels and improves the final quality of the spectra. 
This is especially important for the LH module since it 
contains a large number of such pixels. 

We used Basic Calibrated Data (BCD), processed by 
Spitzer pipeline version S15.3 (low resolution modules: 
SL and LL) and version SI 7. 2 (high resolution mod- 
ules: SH and LH) to build data cubes with CUBISM 
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(|Smith et al.ll2QQ7aD . CUBISM is a tool that combines 
slit observations to create spectral cubes. It also per- 
forms background subtraction, flux calibration and esti- 
mates the statistical uncertainty at each spectral wave- 
length. 

3. ANALYSIS 

3.1. SL Spectral Maps 

Using the SL cubes generated by CUBISM we con- 
structed maps of the 5.5 fivn continuum, and the most 
prominent aromatic features (hereafter PAH) in this 
module (6.2 /im PAH, 7.7 /im PAH and 11.3 /im PAH). 
To build them we used our own IDL procedures. First 
we extracted 1-D spectra from the cubes by averaging 
the spectra over a 2x2 pixel aperture and then by mov- 
ing with a 1 pixel step until the whole field of view was 
covered. The 2x2 pixel extraction box allows us to in- 
crease the S/N of the 1-D spectra, with only a slight 
loss of angular resolution (see below and Appendix [Aj . 
Once we obtained the 1-D spectra, for each pixel we did 
the following. To create the 5.5 /im continuum map, we 
integrated the flux between the rest frame 5.3 /im and 
5.7 /im wavelengths (in this section all wavelengths refer 
to rest frame wavelengths). For the PAH maps, first we 
fitted the local continuum using a linear fit. The con- 
tinuum pivots were set by averaging the spectrum over a 
~0.2 /im interval centered at 5.75 /im and 6.7 /im for the 
6.2 /im PAH, at 7.2 /im and 8.2 /im for the 7.7 /im PAH, 
and at 10.6 /im and 11.8 /im for the 11.3 /im PAH. Af- 
ter subtracting the continuum, we integrated the flux in 
the wavelength ranges between 5.9 and 6.5 /im 7.3 and 
7.9 /im and 10.8 and 11.6 /im for the 6.2, 7.7 and 11.3 
/im PAH features, respectively. We derived the flux un- 
certainty of these maps with the statistical uncertainties 
provided by CUBISM. To avoid noisy pixels, we required 
the feature peak to be at least above la to plot that pixel 
in the map, where a is the continuum noise. The SL 
spectral maps of our sample are shown in Figure [TJ 

The angular resolution (FWHM) of the SL spectral 
maps is ~4.5 arcsec. This is the resolution when we use 
2x2 boxes to construct the maps, as opposed to that 
measured from the standard star observations (see Ap- 
pendix [A]). For the distances of our galaxies corresponds 
to ~1 kpc (the exact scale for each galaxy is given in 
Table [1]) . In Appendix [A] we give full details on how the 
angular resolution of the spectral maps was derived. 

We constructed maps of different PAH ratios by di- 
viding the PAH maps obtained as described above. To 
avoid noisy pixels in the ratio maps we imposed further 
restrictions, and only those pixels with signal-to-noise 
ratio greater than 4 are plotted in the PAH ratio maps. 
Figures[2]and[3]show the 11.3 /im PAH/6.2 /im PAH and 
11.3 /im PAH/7.7 /im PAH ratios, respectively. 

3.2. SH Spectral Maps 

We used the SH cubes to measure the bright- 
est fine structure emission lines: [S IV] 10.51 /im, 
[Ne II] 12.81 /im, [Ne III] 15.56 /im, [S III] 18.71 /im and 
the molecular hydrogen lines, H2 S(2) at 12.3 /im and 
H2 S(l) at 17.0 /im. First we extracted 1-D spectra using 
the same method as for the low r esolution cubes. Then 
we used the IDL package MPFIT ( Mark war dt|l2009D to fit 
simultaneously a first-order polynomial to the continuum 



and a Gaussian to the line. The fit was over a wavelength 
range of ~0.4 /im centered at the line wavelength. The 
spectral resolution, R~600, is sufficiently high to allow 
us to separate the emission lines from the broad PAH 
bands. However, to improve the continuum determina- 
tion, we used extra Gaussian profiles to fit the PAH fea- 
tures close to the emission lines (e.g., the 12.7 /im PAH 
feature next to the [Nell]). To avoid noisy pixels in the 
maps, we required: (1) a signal-to-noise ratio > 2 in the 
continuum range; and (2) the line peak to be at least la 
above the continuum noise for the pixel to be shown in 
the corresponding map. 

We also constructed an 11.3 /im PAH feature map by 
the same procedure as for SL cubes, but with the high 
resolution spectra. A 15 /im continuum was obtained by 
integrating the flux between 14.9 /im and 15.1 /im. We 
required a signal above la to plot a pixel in these two 
maps. 

The SH spectral maps are shown in Figure |4] to- 
gether with the HST/NICMOS continu um-subtracted 
Paa maps of lAlonso-Herrero et al.l (|2006f ). 

We constructed [S III] 18.71 /im/[Ne II] 12.81 /im, 
[Ne III]15.56 /im/[Ne II]12.81 /im and [Ne II]12.81 /im/ 
11.3 /im PAH ratio maps by dividing the flux maps. The 
variation of the PSF size with wavelength within each 
module is negligible (see Appendix [A]), and therefore 
we divided the flux maps directly without attempting 
to correct for this small effect (we only built ratio maps 
of lines observed with the SH module). We imposed 
stricter conditions to plot a pixel in a ratio map than 
those used for the line maps; only pixels with a flux 
measurement > 3a, for both features, were included 
in the ratio maps. The maps of the SH line ratios are 
shown in Figure [5l The angular resolution of the SH 
maps is ~6 arcsec. This corresponds to physical scales 
between ~1 kpc, for the nearest galaxies, and r^2 kpc 
for the most distant ones. 

3.3. Nuclear, Integrated and Extranuclear 1-D Spectra 

We extracted 1-D spectra of the nuclei as well as of the 
integrated galaxy. To extract the nuclear spectra we used 
a 5.5" x5.5" square aperture centered at the nuclear coor- 
dinates (Table [1]). This aperture corresponds to physical 
scales between 1 and 2 kpc, depending on the galaxy. We 
set the orientation of the aperture aligned with the R.A. 
and Dec axis. We used the same aperture for both low 
and high resolution data cubes. To obtain the integrated 
spectrum we summed the flux of the cube. We excluded 
those pixels with the 11.3 /im PAH flux (15 /im contin- 
uum) below la to avoid adding too much noise to the 
low resolution (high resolution) integrated spectra. This 
threshold is the same as that required to plot a pixel in 
the 11.3 /im PAH and 15 /im continuum maps of Figure 
[U and Figure |4] respectively. The physical scale of the 
integration area varies from galaxy to galaxy. The ap- 
proximated physical sizes are listed in Tables [2] and [3] 
for the low and high resolution maps respectively. They 
are of the order of 10x8 kpc. Although the extraction 
regions are slightly different for each module, the flux 
differences are only ~ 10-20%. 

We plot the integrated spectra in Figure [6l The nuclear 
spectra will be presented in Pereira-Santaella et al. (in 
preparation), together wit h the nuclear spectra of the 
rest of the .Alonso-Herrero et al.l (|2Q06.) sample. 
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In Figure [3 (right) we show a comparison of the nu- 
clear and extranuclear low resolution spectra normalized 
at 14 /im for 4 of the galaxies. We selected these galaxies 
because in their spectral maps (left panel of Figure[7]) it is 
possible to identify extranuclear H ll regions and extract 
their spectra. All the spectra, nuclear and extranuclear, 
show bright PAH features and fine structure lines. The 
variations of the spectral features from galaxy to galaxy 
and within a galaxy are relatively small, at least when 
comp ared to those observed in ULIRGs (Armus et al. 
I2QQ7!) or those observed in star-forming galaxies with 
different metallicities (lEngelbracht et al.H2QQ8f ). That is, 
the comparison of the spectra shows that these LIRGs 
have a more or less uniform shallow silicate absorption 
and relatively high metallicity. 

The PAH ratios, PAH equivalent width and silicate 
depth measured in the nuclear and integrated low reso- 
lution spectra are listed in Table O Table [3] and Table |4] 
show respectively the integrated and nuclear fine struc- 
ture line fluxes in the high resolution spectra. 

3.4. Silicate Feature Maps 

Amorphous silicate grains present a broad feature cen- 
tered at 9.7 jiia that can be measured from the SL 
cu bes. We used a m ethod similar to that proposed 
by iSpoon et al.l (|2QQ7f ) to measure the observed silicate 
strength. First, we extracted 1-D spectra from the cube 
as described in Section 13. li For each spectrum, we es- 
timated the continuum expected at 10 /im by fitting a 
power law through the feature-free continuum pivots at 
5.5 /im and 13 /im. Then we compared the estimated 
continuum and the measured flux using the following ex- 
pression: 



5'si In 



fobsjlO /im) 

fcont{lO /im) 
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where ^si is the silicate strength. When the silicate 
strength is negative, the silicate feature is in absorption, 
whereas a positive silicate strength indicates that the fea- 
ture is seen in emission. 

We estimated the flux at 10 /im instead of 9.7 /im (the 
maximum of the feature) to avoid a possible contamina- 
tion by the molecular hydrogen emission line H2 S(3) at 
9.67 /im. 

There are several factors that affect the accuracy of 
the silicate strength measurement. The most impor- 
tant source of uncertainty is the continuum estimate, 
because of the presence of broad PAH bands in the spec- 
tr a. We selected conti nuum pivots similar to those used 
bv ISpoon et al.l ()2007l ). so that we can make meaningful 
comparisons with their results. The Ssi is equivalent to 
the optical depth of the silicate absorption (rsi) defined 
in IRieke & Lebofsky (1985). Alternatively, a more accu- 
rate way to measure the sihcate optical depth is to fit all 
the features (PAHs, emission lines and dust continuum) 
presen t in the spectrum as PAHFIT does (jSmith et al.l 
l2007bl ). However, to avoid ambiguities between the PAH 
strength and the silicate absorption, the entire low res- 
olution spectral range (5.2 /im -38 /im) is desirable. At 
longer wavelengths (>14.5 /im) the pixel size is about 2.5 
times larger than at shorter wavelengths and since we are 
interested in the spatial distribution of the silicate fea- 
ture we decided to use the first method to estimate the 
continuum. To compared these methods, first we fitted 



the integrated spectra using PAHFIT, assuming a fore- 
ground screen dust geometry to obtain the rg.r, and then 
we compared this value with the measured ^si- We found 
that both quantities are well correlated. 

Keeping in mind that the main source of uncertainty 
of the silicate strength measurements is the continuum 
determination, we estimated that the typical values of 
the statistical uncertainty are ~0.05 for the nuclei and 
less than ~0.15 for the extranuclear regions. The maps 
of the sihcate feature strength are presented in Figure [51 

3.5. SH Velocity Fields 

The SH spectra were also used for rough determination 
of velocities. Figure [9] shows an example of line fitting 
for selected regions of IC 5179 to illustrate the detection 
of variations in the line position, even with the relatively 
low spectral resolution of the SH module. 

Since the spectral resolution, R~600, is lower than that 
of most kinematic studies, we tested the validity of the 
velocities with synthetic spectra. We created these spec- 
tra to simulate an unresolved spectral line as seen by the 
SH module and then we tried to recover the original line 
parameters. The equivalent width of the simulated lines 
was similar to the observed ones. We added noise to the 
spectra to match the continuum signal-to-noise ratio of 
the maps. The continuum signal-to- noise ratio, for most 
of the galaxies, is in the range ~60-100 at the nuclei and 
'^3-10 at the external regions. To allow for the effect 
of the telescope pointing uncertainties, we distorted the 
wavelength scale by ^5% (estimated from the pointing 
information in the BCD file headers). We found that 
this effect limits the accuracy of the velocity determina- 
tion to ± 10 km s~^, independent of the signal-to-noise 
ratio of the spectra. For extended sources, which fill uni- 
formly the slit, the wavelengths are more stable, and then 
the uncertainty due to the pointing inaccuracy would be 
lower. The galaxies studied here are not point sources, 
neither uniformly extended sources, thus this uncertainty 
estimation is an upper limit. Including the uncertainty in 
the absolute wavelength calibration (equivalent to 10% 
of a pixel), we found the minimum uncertainty to be ± 
20 km s~^. We therefore conclude that the uncertainty 
in the velocity field is in the range of 10 - 30 km s~^, 
depending on the pointing accuracy achieved and the 
signal-to-noise ratio of an individual spectrum. Thus, 
we consider that variations in the velocity maps larger 
than ~20 km s~^ are likely to be real. Figure [TOl shows 
the velocity fields derived for the [Nell] 12. 81 /im and 
H2 S(l) emission lines; they range up to a total gradient 
of ~200 km s~^ and hence include a substantial amount 
of information about the sources. 

4. THE SILICATE FEATURE 

4.1. Spatially Resolved Measurements 

LIRGs are known to contain highly obscured re- 
gions (Ay ^ 4 — 50 ma^, see e.g., IVeilleux et al.l 119951 : 
Genzel et al. 1998; Alons o-Herrero et al.l I2000L l200i ' 



usually coincident with the nuclei of the galaxies. As- 
suming an extinction law and a dust distribution geom- 
etry the silicate strength can be converted into a visual 
extinction. 

The maps of the silicate feature of our sample are 
shown in Figure [HI In a large fraction of the systems 
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Fig. 6. — Rest-frame SL (left panel) and SH (right panel) integrated spectra of the sample. The Arp 299 spectra was shown in 
lAlonsQ-Herrero et al.l (|20Q9 b') . The fluxes have been scaled for better comparison. The scale factors are, starting from the top galaxy, 
2x10^, 6x10^, 1.2x10^, 1.5x10^, 2x10"^, 3x10^, 7x10^, IxlO^, 
2x10^, 2x10"^, 2x10^, 5x10^, 2x10^, 2, 2xl0-\ SxlQ-^^ 2x10" 
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5x10-^, 2x10-^ for the left panel and IxlO'^, 
lxlO~^ for the right panel. 



(80%) the most obscured regions appear to be coinci- 
dent with the nuclei. In other cases, the nuclear regions 
appear to be slightly less obscured than the surround- 
ing regions (e.g., NGC 7130) or they show a compli- 
cated morphology (e.g., IRAS 17138-1017, NGC 7771, 
IC 4687). We note that the spectra used for measuring 
the strength of the silicate feature are averaged over ~1 
kpc scales, and thus some of the real variations seen on 
smaller scales using ground-based mid-IR observations 
are smoothed out (e.g., the hig hly obsc ured southern nu- 
cleus of NGC 3256, see Diaz- Santos et al. 2008). 

Figure [11] shows the distribution of the silicate 
strengths measured from the spatially resolved maps for 
all the galaxies. They are in the range ~0 to ~ —2.0, 



with a median of —0.61. The strengths of the majority 
of the LIRGs are moderate (6'si ~ —0.4 to — 0.9), and 
intermediate between those observed in starburst galax- 
ies (Brandl et al. 2006) and ULIRGs (Spoon et al. 2(303; 
ISirockv et al.ll2008i) . 

If we assume a foreground screen of dust, the sili- 
cate absorption strength can be expressed in terms of 
the optical extincti on (Ay) using the conv ersion fac- 
tor Ay/^si = 16.6 (jRiel^e fc Lebofskvl [19851 ) . We then 
find that the median Ay of the spatially resolved mea- 
surements is ~10 mag. The most extincted region in 
our sample of LIRGs is the nuclea r region of IC 694, 
also kn own as Arp 299- A (see also >Alonso-Herrero et al.l 
l2009bl ). with Ay -37 mag. The foreground screen as- 
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Fig. 7. — Left: Spitzer/IRS SL spectral maps of the 7.7 /im PAH feature. The image orientation is indicated on the maps for each 
galaxy. The scale represents 1 kpc. The black squares are the positions and sizes of the extraction apertures. Right: Low resolution spectra 
normalized at 14 /im. 
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Fig. 9. The left panel is the map of the [S III] 18.71 /im hne 
flux of IC 5179. The right panels show the flt to this line for three 
selected regions of this galaxy. 



80 




0.0 -0.5 -1.0 -1.5 -2.0 

Ssi 

Fig. 11. — Histogram of the spatially resolved measurements of 
the silicate strength for all galaxies measured in boxes of 2x2 pixel. 
The solid grey line indicates the median silicate strength. The 
dashed grey lines deflne the range including 68% of the data points. 

sumption gives the minimum extinction compatible with 
the measurements. Most of the spatially resolved mea- 
surements (68%) are in the range of -0.4 to -0.9. This 
range is much narrower than that obs erved in the n uclear 
regions of ULIRGs (0 > 6'si > -4) (jSpoon et all [2007). 
The similarity of the silicate absorption strengths among 
our sample members suggests that they may be in part 
a product of radiative transfer in dusty regions. We can 
approximate this possibility by assuming that the dust 
and emission sources are mixed. The resulting median 
Ay is then ~20 magnitudes, about 2 times greater than 
for the foreground screen geometry. 

The nuclear and integrated values of the silicate 
strength are listed in Table El For all the galaxies, we 
find that the integrated strength of the silicate feature is 
smaller than, although comparable to, the nuclear one. 
The nuclear silicate strengths range from -0.4 to -1.9, al- 
though most of them are larger than -1.0. The integrated 
strength is —0.3 < ^si < —1.2, somewhat shallower than 
the nuclear strength. The LIRGs with the largest varia- 
tions between the nuclear and integrated silicate features 
are NGC 3256, Arp 299, NGC 7130, and NGC 7591. 

4.2. The Ssi vs. 6.2 jim PAH EW Diagram 

iSpoon et al.l (2007) presented a diagram using the 
equivalent width (EW) of 6.2 jiia PAH feature and the 
silicate strength to provide a general classification of in- 
frared galaxies. ULIRGs appear distributed along two 
branches. The horizontal branch is populated by unob- 
scured AGNs (QSO and Syl) in the left-hand side of 
the diagram (class lA) and by starburst galaxies in the 



right-hand side (class IC) of the diagram. The diagonal 
branch goes from deeply obscured nuclei (class 3A) to 
unobscured starbursts (class lA). 

In Figure [12] we present a similar diagram for our 
sample of LIRGs. The left panel shows the spatially 
resolved measurements of our sample of LIRGs. Most 
of them appear in this diagram at the region where 
both branches intercept. The nuclear and surrounding 
regions of galaxies classified as Seyferts using optical 
spectroscopy (NGC 5135: Sy2, NGC 7130: LINER/Sy 
and Arp 299-Bl: Sy2, see Table [1]) are located in the 
IB region (intermediate between AGN and starbursts), 
in agreement with their well-known composite nature. 
NGC 6701 and NGC 7591 are classified as composite ob- 
jects which is likely to be a combination of an AGN an d 
star- formation activity (see lAlonso- Herrero et aIll2009aD . 
Both nuclei appear quite close to region IB. The two 
most obscured nuclei in the sample, Arp 299-A and the 
southern nucleus of NGC 3256, and their neighboring 
regions, are in the upper part of the diagonal branch 
populated by ULIRGs, although they do not reach the 
silicate strengths of the deeply embedded ULIRG nuclei. 

The comparison between the nuclear and integrated 
values (right panel of Figure [T2]) shows that the inte- 
grated spectra tend to have lower silicate strengths and 
larger 6.2 jim PAH EWs. We find that on average the 
integrated 6.2 jiia PAH EW is ~30% larger than the nu- 
clear values, whereas the integrated silicate strength is 
~15% lower than that found in the nuclei. That is, the 
integrated values move the galaxies to the pure star- 
burst class region (IC) in this diagram. In the case 
of all the active galaxies of the sample, the nucleus 
is not sufficiently bright as to dominate the integrated 
mid-infrared spectrum and the dominant contribution 
fro m star formation, which is ex tended over several kpc 
(see lAlonso- Herrero et al.1l2006l [2Q0 9b.a) makes the in- 
tegrated spectrum look more starburst-like. 

5. FINE STRUCTURE EMISSION LINES 
5.1. Morphology 

Although the limited spatial resolution (~ 1 kpc) of 
the SH maps does not allow us to study the galaxy mor- 
phologies in great detail, it is possible to analyze general 
trends. The exact spatial distribution of the emission 
lines depends on the physical conditions and the age of 
the dominating stellar population in each region. In Ap- 
pendix [B] we discuss the galaxies individually. 

Figure S] shows the SH spectral maps of the most 
prominent features, as well as those of the 15 jiia con- 
tinuum. The [Nell]12.81 /im, [Ne III] 15.56 /im and 
[S III] 18.71 /im fine structure lines have an overall mor- 
phology similar to that of the 15 /im continuum. We 
detect [S IV] 10.51 /im emission in only five systems. 
In two of them (NGC 7130 and NGC 5135), the 
[S IV] 1 0.51 /im emission mainl y comes from the active 
nuclei (iDiaz-Santos et al1l2010[ ). Our data do not allow 
us to determine whether there is a low surface brightness 
[S IV] 10.51 /im emission, associated with Hll regions, in 
these two galaxies due to the relatively low fiux of this 
line. In the other three cases, NG C 3256, IC 4687, and 
Arp 2 99 (for the last galaxy, see lAlonso-Herrero et al.l 
'2009b), the emission is more extended and we attribute 
it to the Hll regions seen in the Paa and/or Ha images. 
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Fig. 12. — S'si versus the 6.2 /xm PAH equivalent width (EW) diagram, similar to that presented by Spoon et alj (|20Q7I V Left: Spatially 
resolved measurements for our sample of LIRGs using a 2x2 pixel box (3.7 "x3.7 " aperture). The black crosses are the two nuclei of 
Arp 299 (IC 694 and NGC 3690) and the southern nucleus of NGC 3256. Right: Nuclear (5.5 "x5.5 " aperture) and integrated values of 
each galaxy. The small symbols correspond to the nuclear value and the large symbols to the integrated galaxy. Nuclear and integrated 
are connected by a line to guide the eye. In both panels, following Spoon et al. (2007), the parameter space is divided into 6 classes out 
of their 9 classes. Class lA is populated by unobscured AGNs. They have a featureless spectra with a shallow silicate feature. Class IC 
corresponds to pure starbursts (PAH dominated spectra). The spectra of a galaxy in the class IB is intermediate between the AGN spectra 
and the pure starburst. Classes 2B and 2C have a deeper silicate absorption, so do classes 3A and 3B (not plotted here). The majority of 
ULIRGs not dominated by an AGN lie in these classes. 



We find that the neon emission ([Ne II] 12.81 /im and 
[Ne III] 15.56 /im) and the 15 /im continuum are spatially 
resolved. We calculated the ratio between the nuclear 
(circular aperture of radius 2 kpc; depending on the dis- 
tance the aperture radius varies between 1.5 and 3 pix- 
els )EI and the integrated neon line emission and 15 /im 
continuum emission for our galaxies, we also used the 
uncertainty maps to estimate the statistical uncertainty 
of these ratios (Table [5]). On average the uncertainty of 
the ratio is ~6%. These are upper limits to the relative 
contribution of the nucleus to the total emission of the 
galaxy as in some cases the IRS maps do not cover the 
full extent of the galaxy. In general the emission arising 
from the central 2 kpc accounts for less than 50% of the 
total emission for the neon lines and the continuum, in- 
dicating that the star formation is extended over several 
kpc. We note however that the LIRGs in this work were 
selected because they show extended Paa emission (Sec- 
tion [211]), and thus, this result needs to be confirmed for 
a complete sample of LIRGs. 

The combined luminosities of the [Ne II] 12.81 /im and 
[Ne III] 15.56 /im lines are good tracers of the star forma- 
tion rate in galaxies (Ho & Keto 2007), so a good mor- 
phological correspondence is expected between them and 
Ha. To further explore the behavior of the neon emis- 
sion. Figure [13] con ipares the IRS neon spectral maps 
with Ha imagefl of ' Hattori etlo] (|200l ) for NGC 3110 
and NGC 7771. The SH and Ra maps of NGC 3110 show 
emission from the nuclear region and from H II regions in 
the spiral arms. The [Ne III] 15.56 /im emission (also the 
[S III] 18.71 /im emission) from these Hll regions is com- 
parable to the emission arising from the nucleus, while 
[Ne II] 12.81 /im is more concentrated around the nucleus. 

^ The physical size of the nuclear aperture was chosen so it cor- 
responds to the minimum spatial resolution of the most distant 
galaxy in the sample. 

^ The Hq; images were degraded to the Spitzer /IRS SH angular 
resolution, pixel size, and orientation. 



We note, however, that the Ha image is not corrected for 
extinction. This is especially important in the nuclear re- 
gion (^si ~ —0.48, see Table [2] Ay ^ 8 mag) as it is the 
most extincted region in this galaxy. For NGC 7771, 
the [Nell] 12. 81 /im [Ne III] 15. 56 /im and Ha emissions 
show similar morphologies, both in the nuclear ring of 
star formation (seen in the Paa image. Figure [U but 
not resolved by the IRS spectral maps), and the extra- 
nuclear Hll regions, as well as in the diffuse regions be- 
tween them. These two examples illustrate that at least 
from a qualitative point of view and on the scales probed 
here, the Ha, [Ne II] 12.81 /im and [Ne III] 15.56 /im lines 
are tracing the same young ionizing stellar populations 
fsee lHo fc Keto|[2007l for a quantitative assessment). 

5.2. Spatially Resolved Line Ratios 

5.2.1. The [Ne III]15.56 iim/[Ne Il]12.81 iim Line Ratios 

The spectral maps of the [Ne III] 15.56 /im/ 
[Ne II] 12.81 /im line ratios for our sample of LIRGs are 
shown in Figure [5] In most cases the nuclei have the 
lowest [Ne III]15.56 /im/[Ne II]12.81 /im ratios, while 
the circumnuclear regions and H ll regions tend to show 
higher ratios (see also Figure [H]). With the exception 
of three cases with an active nucleus (NGC 3690, 
NGC 5135 and NGC 7130), the nuclear values of the 
[Ne III]15.56 /im/[Ne II]12.81 /im ratio are ah in the 
range of 0.06 - 0.10 (with three exceptions near 0.2, 
IC 694, IRAS17138-1017 and IC 4687 - see Table p. 
Excluding the galaxies hosting an AGN, the median 
nuclear ratio is 0.08. The extra- nuclear values are 
between 0.06 and 1, while the integrated values (see 
Table [3]) are between 0.08 and 0.5. T hornlev et al.l 
(|2000[ ) find a similar range for the integrated ratios 
for an additional four LIRGs. Thi s ratio range s froni 
0.02 to 6 in starbu rst galaxies (|Verma et al.l [20031 : 
[Dale et al.[[2QQ6[ . [2009[ ). For the Hll galaxies, the nuclear 
[NeIIl]15.56 /im/[Ne II]12.81 /im ratio is a factor -2-3 
smaller than those measured in the extranuclear regions. 
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Fig. 13. Comparison of the [Ne Il]12.81 /im and [Ne III]15.56 /im maps with Ha images from lHattori et al.l pOOl ) for NGC 3110 and 
NGC 7771. The image orientation for each galaxy is shown on the Ho; image. All the maps are shown in a square root scale. 



A similar behavior of increasing [Ne III] 15.56 jam/ 
[Ne II] 12.81 /im ratios with increasing galactocentric dis- 
tances has been obse rved in the nearby galaxies NGC 253 
(|Devost et al.l I2QQ41 ) and M82 (Beirao et alj |2008i), al- 
though the physical scales probed are different (^80 pc 
in these two galaxies, ~1 kpc in our sample of LIRGs). 
Galactic H ll regions follow the same trend, with H ll re- 
gions at increasing distances from the galactic nucleus 
showing larger [Ne III] 15.56 /im/[Ne II] 12.81 /am ratios 
(jGiveon et al. 2002). 

Note that we did not correct this ratio for extinction. 
This correction however, only produces variations in the 
ratio of 5 — 15% for the typical extinction of our sample 
{Ssi ^ —0.6, that is Ay ~ 10 mag). The correction 
(positive or negative) depends on the adopted extinction 
law (see Farrah e t al.ll2007[ ). 

The strength of [Ne III] 15.56 /im relative to 
[Ne II] 12.81 /im decreases with increasing metallic- 
ity, increasing age of the ste har populatio n and 
increasing nebular density (e.g., iSnijders et al.l [2007). 
In accordance with the luminosity- and mass-metallicity 
relations, LIRGs are characterized by uniformly rel- 
atively high metallicity (in many cases super-solar), 
with some indications of sli ghtly (factor of tw o) re- 
duced levels i n their nuclei (|Rupke et al.l 120081 ) . The 
ISnijders et al.l ([2007) models predict a ~3 times larger 
[Ne III] 15. 56 /im/[Ne II] 12.81 /im ratio for a given age 
and density as the metallicity decreases from super-solar 
(2 Zq) to solar. Thus, metallicity-related effects are 
unlikely to be the dominant cause of the relative increase 
of [Ne III] 15.56 /im strength outside the central regions 
of these galaxies. They may in fact to be overcome by 
some stronger mechanism. 



The trend with stellar population age is strong. The 
expected [Ne III] 15.56 /im/[Ne II] 12.81 /im ratio of an in- 
stantaneous burst of star forma tion falls below 10~^ 
for ages >6 Myr ( Rigbv fc Riekd 120041 ) . iThornlev et all 
(2000) showed that the range of the [Ne III] 15. 56 /im/ 
[Ne II] 12.81 /im ratio in starburst galaxies is compatible 
with models for bursts of star formation. Howe ver, from 
their models and those of iRigby fc Riekd (|2004l ). the nu- 
clear values for LIRGs are below the range that would be 
expected from such star-forming episodes. The observed 
nuclear ratios could only be explained in terms of the 
stellar population if: (1 ) the IMF is tr uncated at ~30 
M0 in the LIRG nuclei (jThornlev et all 2000, Figure 6); 
or (2) the star formation rate in all the LIRG nuclei has 
decreased rapidly over the past 10-20 Myr (which seems 
contrived) . 

Another explanation is that due the high densities in 
the nuclear regions the [Ne III] 15.56 /im/[Ne II] 12.81 jam 
ratio is suppressed, including the possibility that a 
fraction of the most massive st ars are hid den in 
ultra-compact Hll regions (Rigby fc Rieke'^^OOi). The 
[Ne III] 15.56 /im/[Ne II] 12.81 jam ratio decreases by ap- 
proximately a fact or of two if the dens ity increases from 
10^ to 10^ cm-^ (jSniiders et al.ll2007[ ). We explore the 
density effects on the neon line ratio in the following sec- 
tion. 

Before doing so, we point out that a few galaxies 
(NGC 3690, NGC 5135 and NGC 7130) have an opposite 
trend in that the [Ne III]15.56 /im/[Ne II]12.81 /im ratio 
becomes larger in their nuclei. We attribute this behav- 
ior to the presence of active nuclei. These three galaxies 
are the only ones classified as Seyfert type and it is likely 
that the extra hard-ionizing radiation from the AGN is 
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Fig. 14. — Range of the observed [Ne III] 15.56 /xm/ 
[Nell] 12. 81 fim ratio fr om the spectral maps. The galaxy 
symbol (as in Figure I12|) indicates the value of the ratio in the 
nucleus. The horizontal mark is the median of the ratio of each 
galaxy. The solid gray line is the the median of the medians and 
the dotted line is the median of the nuclear ratios. The black line 
indicates the typical uncertainty of the ratios. 

responsible for the high, three times larger than the me- 
dian, [Ne III] 15.56 /im/[Ne II] 12.81 /im nuclear ratio. 

5.2.2. The [S III]18.71 ^im/[Ne Il]12.81 /am Line Ratios 

The excitation potentials of the [Ne II] 12.81 /im 
and [S III] 18.71 /im emission lines are 21.6 eV and 
23.3 eV, respectively. Therefore, the [S III] 18.71 /im/ 
[Nell] 12. 81 /im ratio is almost insensitive to the hard- 
ness of the radiation field. On the other hand, this ra- 
tio is a good tracer of density in the range 10^ — 10^ 
cm-^, with lower [S III]18.71 /im/[Ne II]12.81 /im ratios 
indicating larger electron densities (Sniiders et al. 2007). 

Our spatially resolved measurements show that this 
ratio ranges from ~0.1 to ~1.25 (Figure [15]). As- 
suming a 5 IVIyr old stellar population, an interme- 
diate ionization parameter (q = 1.6 x 10^) and so- 
lar metallicity (see caption of Figure [T5]) . these ratios 
corres pond to electron de nsities between 10^ and 10^ 
cm-^ (|Sniiders et al.ll2007D . Like the [Ne III] 15.56 /im/ 
[Ne II] 12.81 /im ratio, for a given galaxy the lowest val- 
ues of the [S III] 18.71 /im/[Ne II] 12.81 jam ratio occur in 
the nuclei (in most cases the nuclear value is between 
'^0.13-0.25), while higher ratios are found around them 
(Figure [5]). This indicates, as expected, that the nuclei 
are denser than the extranuclear Hll regions. According 
to this spatial distribution the integrated ratios (~0.2- 
0.5) are ~30% larger than the nuclear ratios. 



The left panel of Figure [16] shows the [S III 



18.71 /im/ 
15.56 /im/ 



Ne II] 12.81 /im ratio versus the [Ne III 
Nell] 12. 81 /im ratio from our spatially resolved mea- 
surements. As can be seen from this Figure, there is a 
tight correlation between the two line ratios if we exclude 
those galaxies hosting an active nucleus (NGC 3690, 
NGC 5135 and NGC 7130). The best linear fit to the 
data (excluding the AGNs) is: 

log{y) = 0.19 ± 0.02 + (0.81 ± 0.03) x log{x) (2) 

Where y is the [S III] 18.71 /im/[Ne II] 12.81 /im ratio 
and X the [Ne III] 15.56 /im/[Ne II] 12.81 /im ratio. In 
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Fig. 15. Range of the observed [S III] 18.71 /im/[Ne II] 12.81 /im 
ratio from the spectral maps. The symbols are as in Figurefl4l The 
density label corresponds t o the [S III] 18.71 /im/[Ne II] 12.81 /im 
ratio predicted by the Sniid ers et al.l (|2007) models for solar metal- 
licity, q = 1.6 X IQS and age = 5 Myr. The [Ne III] 15. 56 /im/ 
[Ne II] 12.81 /im ratio predicted using these parameters is ~0. 1-0.2 
which is in agreement with the observed ratio. 

the right panel of Figure [16] we plot the nuclear ra- 
tios of our LIRGs together with those of extranuc lear 
Hll regions and nuclei of nearby galaxies (Dale et al.l 
120091 ), S eyfert galaxie s (iTommasin et all |2008[ ) and 
quasars (|Veilleux et al.l l2QQ9f ). Both dia srams show 
that, in general, for a given density (as traced by the 
[S III]18.71 /im/[Ne II]12.81 /im ratio) the AGNs show 
larger [Ne III] 15.56 /im/[Ne II] 12.81 /im ratios than star- 
burst galaxies and Hll regions. The most straight- 
forward explanation is that the ionizing photons from 
the active nucleus increase the radiation hardness (the 
[Ne III] 15.56 /im/[Ne II] 12.81 /im ratio). However some 
of the Seyfert galaxies and quasars are in the diagram 
close to Hll regions. These may be galaxies harboring 
an AGN as well as star formation. IVloreover the posi- 
tion of the AGNs a nd starbursts i n Fi gure [T6l is consistent 
with models (^Groves et all 120041 : IPopita et al...2006. ). 

5.2.3. The [S IV]10.51 /im to [Ne III]15.56 /im and to 
[S III]18.71 iim Line Ratios 

The [S IV] 10.51 /im and [Ne III] 15.56 /im lines arise 
from ions with a similar excitation potential, 35 eV and 
41 eV, respectively. While the [Ne III] 15. 56 /im line is 
detected in all galaxies, the [S IV] 10.51 /im line is de- 
tected in just six galaxies (40% of the sample). If the 
[S IV] 10.51 /im emission arises from silicate dust embed- 
ded regions, it can be heavily affected by extinction as 
it is inside the broad 9.7 /i m silicate feature. This seems 
to be the case in ULIRGs ( Farrah et aL|[2QQ7 ). However, 
the depth of the silicate features in our sample of LIRGs 
are moderate, and generally much lower than in ULIRGs 
probably indicating lower extinction (see Section |4]). 

Usin g the photoionization models of I Sniiders et al.l 
(|2007l ), we can infer the expected range of 
[SIV]10.51 /im/[NeIIl]15.56 /im ratios. The values 
for the range of densities, ionization parameters and 
stellar population ages covered by the models, excluding 
the most extreme cases, are compatible with the upper 
limits we measure. Moreover, we do not find any 
correlation between the [S IV] 10.51 /im/[Ne III] 15.56 /im 
ratio and the observed silicate strength (^si). Thus it is 
likely that the non-detection of the [S IV] 10.51 /im line 
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Fig. 16. — Left panel: [S III] 18.71 /im/[Ne II] 12.81 /im ratio vs. [Ne III] 15.56 /im/[Ne II] 12.81 /im ratio from the spatially resolved maps 
(Each point corresponds to a resolution element including the nuclei). Right panel: Nuclear ratios of our sample of LIRGs. We can compare 
them with: Hll regions (grey stars) and star-forming galaxies (red circles) of Dale et al. (2009) , Seyfert galaxies of Tommasin et al. (2008) 
(green circles) and quasars of Veille ux et al.. (^200^) (blue diamonds). The LIRG symbols are as in Figure fl4] The black dashed line is the 
best fit to the spatially resolved data of the LIRGs excluding those galaxies containing an active nucleus which lie below the correlation. 
We did not apply an extinction correction. 
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Fig. 17. [S IV] 10.51 /im/[S III] 18.71 /im ratio vs. 
[Ne III] 15.56 /im/[Ne II] 12.81 /im ratio from the spatially re- 
solved maps. We did not apply an extinction correction. The 
galaxy symbols are as in Figure 1141 When [S IV] 10.51 /xm 
is not detected we represent the upper limit of the nuclear 

[S IV] 10.51 i2m/\S III] 18.71 rati o The dashed line is the 

correlation found bv I Gordon et aP (|2008l ) for MlOl Hll regions 
and starburst galaxies. The black arrow indicates the extinction 
correction of Ay = 10 mag, assuming a screen geometry. The 
orange and light green points above the correlation corresponds to 
galaxies harboring an AGN (see Section r5. 2. 3|) . 



is because its flux is under the detection limit. 

The [S IV]10.51 /im/[S III]18.71 /am and the 
[Ne III]15.56 /im/[Ne II]12.81 /im ratios are weh cor- 
related and can be u sed to estiniate th e hardness 
of the radiation field (IVerma et al.l I2QQ3I : [Dale et al.,, 



I2QQ6I : iGordon et aTll2QQ8[ ). Figure [HI shows that there 
is a good correlation between the [S IV] 10.51 jam/ 
[S III]18.71 /im and [Ne III]15.56 /im/[Ne II]12.81 /im 
ratios of the spatially resolved measurements for our 
star- forming LIRGs, and inde ed, they app r oxima tely 
follow the correlation found by IGordon et al.l (|2QQ8[ ) for 
MlOl Hll regions and starburst galaxies. The upper 
limits to the [S IV] 10.51 /im/[S III] 18.71 /im ratio are 
also consistent with this trend. 
Two of the galaxies classified as active, NGC 5135 and 



NGC 7130, lie above the IGordon et aP (l2008l) c orrela- 
tion. A similar trend was found bv IPale (j2006h 
for AGN in the [S IV] 10.51 /im/[S III]33.48 /im versus 
the [NeIIl]15.56 /im/[Ne II] 12.81 /im diagram. A pos- 
sible explanation is that extra [S IV] 10.51 jam emission 
is produced by the active nucleus. High spatial resolu- 
tion ground-based observations of these two LIRGs in- 
dicate that most of the [S IV] 10.51 /im emission comes 
from the nucleus, whereas this Hue is not detected in 
the surrounding star-forming regions of these galaxies 
(|Diaz-Santos et al. 2010). For the other two galaxies 
NGC 3690 (classified as Seyfert), and NGC 7591 (classi- 
fied as composite, intermediate between H ll and LINER) 
the [S IV] 10.51 /im line is more likely to be affected by 
extinction, as both have nuclear Ssi ~ —0.8, and thus the 
observed [S IV] 10.51 /im/[S III] 18.71 /im ratios are lower 
limits. 

6. PAH FEATURES 

Models predict that the relative st rength of the PAH 
bands depends on the ionization state (jPraine fc Lill200l] : 
iGallianol [2006h . The 6.2 /im PAH and 7.7 /im PAH 
bands dominate the emission of ionized PAH, whereas 
the 11.3 /im PAH is more intense for neutral PAHs. 
More complete models of PAH feature b ehavior have 
been published recently bv iGalHano et al.l (2008), who 
explore the effects of the radiation field hardness and 
intensity as well as the size distribution of the PAH car- 



6.1. Morphology 

Figure [T] shows flux maps of three brightest PAH fea- 
tures present in the SL module at 6.2 /im, 7.7 jam and 
11.3 /im, together with the 5.5 /im continuum maps. Ta- 
ble [5] shows the ratio between the flux arising from the 
central 2 kpc and the total emission. The 11.3 /im PAH 
is clearly more extended than the 5.5 /im continuum in 
60% of the galaxies , arising; also f rom more diffuse regions 
(see also Diaz- Santos et al]l201G) . However, the nuclear 
contribution to the integrated 6.2 /im PAH emission is 
comparable to that at the 5.5 jam continuum, and only 
in 2 cases (NGC 7130 and NGC 7591) the 6.2 /im PAH 
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emission is more extended than the continuum. When 
comparing both PAHs, the 11.3 /am PAH emission is 
more extended than the 6.2 /im PAH emission in about 
40% of the galaxies, while the opposite is never the case. 

A direct comparison of the PAH ratios with fine struc- 
ture line ratios is not possible because the low spectral 
resolution (R~60-120) is not sufficient to separate the 
emission lines from the PAHs. Alternatively we used 
high resolution (R^600) data to measure the fine struc- 
ture lines and the 11.3 /im PAH (Figure H]). 

In general, the 11.3 /im PAH emission appears to 
be more extended than the [Ne U] 12.81 /im line. 
The PAH emission is also more extended than the 
[Ne n] 12.81 /im+[Ne HI] 15.56 /im emission which, as dis- 
cussed in Section O traces recent star formation. PAHs 
have also been used to measure the star formation, 
since they are a proxy, although n ot perfect, for the to- 
tal inf rared luminosity (Pe eters et al.ll2QQ3 : iSmith et al.l 
l2QQ7b[ ). From this comparison, it is clear that the 
11.3 /im PAH emission does not trace the same stars as 
the [N e U] 12.81 jam emission line. In fact, iPeeters et al.l 
(|2QQ4 ^ found that the 6.2 /im PAH is not a good tracer 
of massive star formation (O stars) as seems to be the 
case of the 11.3 jam PAH. 

6.2. Extinction Effects on the PAH Ratios 

As discusse d above, models (|Draine 
iGalliand l2QQ6f ) predict that neutral PAHs show larger 
11.3 /im PAH to 6.2 /im PAH and 11.3 /im PAH to 
7.7 /im PAH ratios than ionized PAHs. Figures [2] and 
[3] show the SL spectral maps of the 11.3 /im PAH/ 
6.2 /im PAH and 11.3 /im PAH/7.7 /im PAH ratios, re- 
spectively. The 11.3 /im PAH/6.2 /im PAH ratio spans a 
factor of 2, from -0.45 to -1.15 and the 11.3 /im PAH/ 
7.7 /im PAH ratio range is — 0.3 — 0.9. The se ranges are 
similar to those found in starburst galaxies (iBrandl et al.l 
12006 ). 

The 11.3 /im PAH is affected by the broad 9.7 /im 
silicate feature, thus in highly obscured regions the 
the 11.3 /im PAH/6.2 /im PAH and the 11.3 /im PAH/ 
7.7 /im PAH ratios can be underestimated. An example 
of this is the southern nucleus of NGC 3256, one of the 
most obscured regions (^si = -1.4) in the sample. Its 
11.3 /im PAH/6.2 /im PAH ratio is lower than that in 
the surrounding regions. 

To explore the effects of the 9.7 jam silicate feature 
on the PAH ratios, we compare the 11.3 jam PAH/ 
6.2 /im PAH ratio and ^si for the spatially resolved mea- 
surements of our LIRGs in F igure [T8l f left panel) . A sim- 
ilar figure was presented by IBrandl et al.l (|2006[ ) for the 
nuclear regions of starburst galaxies. They found that 
the extinction can change the relative strength of the 
PAHs by up to a factor of 2. We find, however, that the 
median value of the 11.3 /im PAH/6.2 /im PAH ratio is 
approximately constant for silicate strengths weaker than 
— —1.0, and only deeper silicate absorptions seem to af- 
fect the 11.3 jam PAH/6.2 /im PAH ratio. Some caution 
is needed here because there are very few regions in our 
sample with silicate strengths stronger than —1.0. 

As an alternative approach, we corrected these PAH 
ratios for extinction using the spatially resolved measure- 
ments of the ^si- We assumed a foregrou nd screen dus t 
geometry and the extinction law of Smith et al.l ()2007bl ). 
For a typical value of the silicate strength in our sample 



of Ssi = —0.5, this extinction reduces the 11.3 /im PAH 
to 6.2 /im PAH and to 7.7 /im PAH ratios by - 15%. 
If we increase the silicate strength up to —1.0 the re- 
duction is — 30%. These are much less than the factor 
of —2 variations in the 11.3 /im PAH/6.2 /im PAH ratio 
observed in some galaxies (e.g., NGC 3256, NGC 3110, 
NGC 6701), at an almost constant value of 5'si. Thus, we 
conclude that the observed variations in the PAH ratios 
are real, and do not significantly depend on the extinc- 
tion, at least for regions in our sample of LIRGs with 
relatively shallow silicate absorptions (6'si > — 1)- 

6.3. Nuclear vs. Integrated PAH Ratios 

The right panel of Figure [18] compares the 
11.3 /im PAH/6.2 /im PAH ratio vs Ssi values of the 
nuclear spectra and the integrated spectra. The gen- 
eral trend is that the integrated spectra have shal- 
lower silicate absorptions (—15% lower ^si, on average) 
and slightly larger (— 20%, on average) 11.3 /im PAH/ 
6.2 /im PAH ratios. 

In Figure [19] we represent the 11.3 /im PAH/ 
6.2 /im PAH ratio versus 7.7 /im PAH/6.2 jam PAH ra- 
tio, comparing the nuclear and the integrated values of 
each galaxy. The distribution of the values of the ra- 
tio is approximately uniform. However it is interest- 
ing to note two clearly different behaviors. The ma- 
jority of the galaxies show a small variation of the 
7.7 /im PAH/6.2 /im PAH ratio and a somewhat larger 
11.3 /im PAH/6.2 /im PAH ratio in the integrated spec- 
trum. This behavior is expected, since the 7.7 /im PAH 
and 6.2 /im PAH are related to ionized PAHs and the 
11.3 jam PAH to neutral PAHs. The ionized PAHs emis- 
sion is concentrated in the nuclear regions, but the 
there is a substantial amount of the 11.3 /im PAH emis- 
sion arising from diffuse regions throughout the galax- 
ies (see Section lO]) . This explains the larger integrated 
11.3 /im PAH/6.2 /im PAH ratios of ah the LIRGs (see 
Table [2]). 

Some galaxies, however, show a completely differ- 
ent behavior, with an almost constant 11.3 jim PAH/ 
6.2 /im PAH ratio but a lower 7.7 jim PAH/6.2 jim PAH 
ratio in the integrated spectrum. The most extreme 
cases of this behavior are NGC 7130 and NGC 7591. 
IGalliano et al.l ((20081) showed that the 7.7 jim PAH/ 
6.2 jim PAH ratio is very sensitive to the size distribu- 
tion of the PAHs, with smaller PAH molecules producing 
lower 7.7 jim PAH/6.2 jim PAH ratios. One possibility 
is that smaller PAH molecules are destroyed in harsher 
environments, such as those of Seyfert and LINER nuclei. 
The destruction of small PAHs might suppress the nu- 
clear 6.2 /im PAH emission. This may lead to the differ- 
ent behavior observed in the AGN/LINER nuclei in Fig- 
ure [19] Moreover, the absolute value of the ratio seems 
to give little information about the nature, starburst or 
AGN, of the nuclei and we can only distinguish them 
by the relative variation between nuclear and integrated 
ratios. 

6.4. The Effects of the Radiation Field 

Since the [NeIIl]15.56 /im/[Ne II] 12.81 jim ratio is a 
proxy for the hardness of the radiation field, we can com- 
pare the spatial variations of this ratio with those of the 
11.3 jim PAH/6.2 jim PAH to ascertain if the radiation 
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Fig. 18. — 11.3 /im PAH/6.2 /im PAH ratio vs. silicate strength (.Ssi). Left: Spatially resolved measurements using a 2x2 pixel box. 
The symbols are as in Figure [M] The solid gray line is the median and the dashed gray lines encircle the region containing 68% of the data 
points. These are calculated in intervals of 0.25 of the Sqi. The top histogram is the distribution of the the 11.3 /im PAH/6.2 /im PAH 
ratio for the whole sample of LIRGs, in gray are marked the median and the 1-cr. Right: Same as left panel, but comparing the nuclear 
(small symbols) and integrated (large symbols) values. For each galaxy the values are connected by a line to guide the eye. 

field has an effect on the PAH ratios. A comparison be- 
tween Figures [2] and [5] indicates that the 11.3 /im PAH/ 
6.2 /im PAH ratio does not decrease necessarily in re- 
gions of enhanced [Ne III] 15.56 /im/[Ne II] 12.81 jam ra- 
tios. 

On the other hand, the morphology of the 
11.3 /im PAH/7.7 /im PAH ratio (Figure 0j resembles 
that of the [Ne III] 15.56 /im/[Ne II] 12.81 /im ratio. How- 
ever this resemblance is the opposite of what one would 
expect. Both ratios present their minima in the nu- 
clear regions for most of the galaxies. Thus while the 
observed nuclear [Ne III] 15.56 /im/[Ne II] 12.81 /im ratio 
would indicate lower radiation hardness in the nuclei, 
the PAH ratio would suggest a larger ionization state. 
Again, it is likely that the typical radiation field hard- 
ness observed in our sample of LIRGs is not sufficiently 
large as to affect the PAH emissi on, as is the case of 
the low metallicity galaxies of Wu et all ([2006). An- 
other possible interpretation is that the PAHs and the 
neon emission are related to different physical processes 
that due to our limited spatial resolution we cannot sep- 
arate. As a result, the bulk of the 7.7 /im PAH emis- 
sion may be produced in star- forming regions, while the 
11.3 /im PAH emission also comes from diffuse regions. 
Thus, low values of the 11.3 /im PAH/7.7 /im PAH ratio 
trace Hh regions, and larger values indicate less recent 
star fo rmation. The l atter interpretation is similar to 
that of iGalliano et all (j2008D . They find that the lowest 
values of the 11.3 jam PAH/ 7. 7 jam PAH ratio are asso- 
ciated with the brightest PAH emission. However, in 
their galaxies, the 11.3 /im PAH/6.2 /im PAH ratio fol- 
lows the same trend, whereas this is not the case for the 
LIRGs studied here. Alternatively, as discussed in the 
previous section, the size of the grains may affect this 
ratio. If the smaller grains are favored in and around the 
nucleus, the 7.7 jam PAH emission would be increased 
respect the 11.3 /im PAH emission. As consequence the 
11.3 /im PAH/ 7. 7 /im PAH ratio would be lower in the 

n ucleus. 

iBrandl et al.l ()2006l ) studied the dependency of the EW 



Fig. 19. 11.3 /im PAH/6.2 /im PAH ratio vs. 7.7 /im PAH/ 
6.2 /im PAH ratio for nuclear and integrated values. Symbols are 
as in the right panel of Figure fTSi 
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of the PAHs with the radiation ha rdness for a sample of 
high metallicity starburst galaxies. IWu et ahl (^006) car- 
ried out a sim ilar study, but for low metallicity galaxies. 
Whereas iBra ndl et al. (2006) found that this dependency 
does not exist, Wu et al. (2006) found the opposite, that 
is, the harsher radiation of low me tallicity galaxi es field 
lowers the PAH EW. Eventually iGordon et al.l (j2008) 
showed that the reason for the discrepancy is the dif- 
ferent range of the radiation field hardne ss in hi^h and 
low metallicity galaxies. In particular, IGordon et al.l 
(^2008) found that for ionizations above a certain thresh- 
old there is a correlation between the radiation hardness 
and the PAH EW, but for lower ionizations this depen- 
dency disappears. We show in Figure [20] the EW of the 
11.3 fim PAH and the [Ne III]15.56 /im/[Ne II]12.81 /am 
ratio of those galaxies not classified as AGN. The ion- 
ization and metallicit y of our LIR Gs are similar to those 
galaxies studied by Brandl et a l.l (2006) and, similarly 
to their findings, there is no correlation between these 
parameters in our sample of LIRGs. 

6.5. Dependency on the Age of the Stellar Population 

Since the [Ne II] 12.81 /im emission traces young stel- 
lar populations (<10 Myr), and since the 11.3 jam PAH 
feature can a lso be excited by older populations (see 
iPeeters et all 2004 and Sect ion [61]), the [Ne II] 12.81 /im/ 
11.3 /im PAH ratio can give us an indication of age 
variations throughout the galaxies (jPfaz- Santos et all 
I2010r ). As can be see from Figure[5l the [Ne II]12.81 /im/ 
11.3 jam PAH ratio ranges from ~0.1 to ^0.6 for our 
sample of LIRGs. In some LIRGs (e.g., NGC 6701, 
NGC 7130, NGC 7771, IC 4687) the [Nell]12.81 /im/ 
11.3 /im PAH ratio appears enhanced in the high sur- 
face brightness Hll regions (see the HST/NICMOS Paa 
images in Figure [4]) with respect to regions of more dif- 
fuse emission. In other galaxies the interpretation might 
not be as straightforward. In some cases the IRS angu- 
lar resolution is not sufficiently high to isolate the cir- 
cumnuclear Hll regions (e.g., NGC 5135), and the age 
differences may be smoothed out. 

7. MOLECULAR HYDROGEN 

We constructed maps of the molecular hydrogen emis- 
sion using the pure rotational lines H2 S(l) at 17.0 /im 
and the H2 S(2) at 12.3 /im (FigureH]). These lines trace 
warm (300 K < T < 1000 K) molecular gas in the in- 
terstellar medium. The H2 S(l) morphology appears, 
generally, different from that of the fine structure lines, 
which are associated with the ionized gas medium. There 
are also differences in the velocity fields that will be dis- 
cussed in Section [8] Using the same procedure as de- 
scribed in Section ISTT] we find that the H2 S(l) emission 
is more extended than the 15 /im continuum in ~60% 
of the sample. In most of these galaxies the size of the 
H2 S(l) emitting region is larger than that of the fine 
structure lines (see Table [5]). 

7.1. Temperatures, Column Densities and Masses 

We extracted the nuclear spectra of all galaxies using 
a 13.4" X 13.4" aperture. This aperture corresponds to 
physical scales between 3 and 5 kpc depending on the 
galaxy. This aperture size is different from that used in 
the other Sections because now we are limited by the 



lowest spatial resolution of the LH module, where the 
H2 S(0) line at 28.2 jam lies. The ffuxes of the H2 lines 
measured from these spectra are listed in Table [6] 

The temperature of the warm m olecular hydrogen was 
calculated as described by Ro ussel et a l. (200 ^. We 
assumed an ortho- to-para ratio = 3 (iRigopoul ou et al.l 
I2OO2I : [Roussel et al. 2007), and thus this assumption may 
introduce some uncertainty as it is used in the tempera- 
ture calculations. The column densities and masses are 
highly dependent on the H2 lines used in these calcula- 
tions. It is preferable to use the lines H2 S(0) and H2 S(l) 
because they come from the lowest rotational levels of the 
hydrogen molecule, and therefore they trace the molecu- 
lar hydrogen gas at the lowest tem peratures which dom - 
inates determination of the masses (lR.oussel et al.ll2007D . 
However the H2 S(0) line is not detected in all galax- 
ies, and to allow meaningful comparisons we estimated 
the column densities and masses with the H2 S(l) and 
H2 S(2) lines. The results are presented in Table [3 The 
masses estimated from the H2 S(l) and H2 S(0) lines are 
in the range 0.4 to 3 xlO^M©. These values ar e com - 
parable to those found in ULIRGs (Higdon et al. 20 Q6D, 
local starbursts and Seyfert galaxies (Rigopoulou e t al.l 
I2OO2I : iTommasin et al.ll2008D . It should be noted that 
these masses correspond to a beam size of 13.4" x 13.4" 
which does not include all the molecular hydrogen emis- 
sion. Then the mass values Hsted in Table [7] should be 
considered as lower limits to the total H2 mass of the 
galaxy. For instance the integrated fluxes of the H2 S(l) 
line (Table [3]) are a factor of between two and three larger 
than those used in the mass calculations (Table [6]). Thus 
assuming that the physical conditions of the H2 in the 
outer parts of the galaxy are similar to those in the inner 
~3-5 kpc, the total mass would be at least a factor of 
~2-3 larger than the values listed in Table [7] 

7.2. Excitation Mechanism and Relation to the PAH 
Emission 

Studies of the molecular hydrogen emission of galac- 
tic PDRs indicate that the form ation of H2 occurs 
in the surfaces of PAHs (see e.g., iHabart et al.l [20031 : 
IVelusamv fc Langer. 20 08) as both the PAH carriers and 
molecular hydrogen can be excited by UV radiation. 
Thus the H2 and the PAH emission morphologies are cor- 
related in these r egions. A similar correl ation is seen in 
external galaxies (Rigo poulou et al.] [20021 : iRoussel et al.l 
2007). However, H2 c an be excited by o ther mechanisms, 
such as X-rays fe. g.. iLepp fc McCrav 1983) and shock 
fronts fe.g.. Hollenbach fc McKee 1989 ). 

To explore the morphological differences we con- 
structed maps of the H2 S(l)/11.3 jam PAH ratio. The 
range of this ratio for each galaxy is shown in Figure 
[2TI This ratio varies by a factor of 10 in our sample of 
LIRGs, although, as can be seen in the Figure 122} 68% of 
the values are in a narrower range, from 0.028 to 0.073. 
Figure [23] shows the ratio maps of those galaxies with 
the most noticeable morphological characteristics. 

The 11.3 /im PAH feature falls close to the the 
9.7 /im silicate absorption and, although the H2 S(l) is 
also close to the sihcate absorption at 18 /im, extinc- 
tion variations could increase artiflcially the H2 S(l)/ 
11.3 /im PA H ratio. We use d the extinction curve 
adopted by [Smith et al.l (|2007bl ) to quantify this effect. 
If we assume a mixed model for the dust distribution, 
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the ratio can be increased by up to a factor of ~1.4, 
for the typical silicate absorption observed in our sam- 
ple of LIRGs. Which is much lower than the observed 
variation. Choosing a screen geometry, to increase the 
ratio by a factor of 2 we would need an optical extinc- 
tion of Ay ~100 mag, which is much larger than our 
estimates (Section H]). Therefore, the effect of the extinc- 
tion is small compared to the observed variations of the 
H2 S(l)/11.3 /im PAH ratio. 

We can try to distinguish among the different excita- 
tion mechanisms by comparing the 11.3 /im PAH mor- 
phology with that of the H2 emission (Figure |4]). It is 
clear that the morphologies are different. The H2 S(l) 
emission is more extended than the PAH emission, al- 
though the maxima of both emissions are located in the 
nucleus for most of the galaxies. The only exception 
is NGC 3256 where the H2 S(l) peaks at the southern 
nucleus while the PAH maximum is located at north- 
ern nucleus. We note, however, that the southern nu- 
cleus of this galaxy is v ery obscured (see Figure [HI and 
iDiaz-Santos et al.ll2QlQf ). and the 11.3 /im PAH map is 
not corrected for extinction. This suggests that the warm 
H2 in LIRGs is produced not only in PDRs, but that the 
other mechanisms also play a role. 

We attribute the increased values of the H2 S(l)/ 



11.3 /im PAH ratio, not to a deficiency of PAH emission, 
but rather to an excess emission by H2 , excited by X-ray 
or shocks. For instance, in the merger system NGC 3256, 
the largest value of the ratio occurs where the compo- 
nents of the merger are colliding. The H2 line ratios 
are inconsistent with X-ray excitation (I Shaw et al. 2005). 
Thus, it is likely that there is some contribution from 
shocks to the H2 , as appears also t o be the case for some 
regions in NGC 3690 and IC 694 (|Alonso-Herrero et al.l 
[2009bh . In many cases, the hypothesis of shock ex- 
citation is suppo rted by large velocity dispersions (see 
IGarcia-Marinet al.l 1200 6 for NGC 3690/IC 694 and the 
following section for additional examples). 

The galaxy IC 5179 is interesting because its veloc- 
ity field (see next section) is highly ordered, so strong 
shocks are unlikely, and its X-ray luminosity is low. Thus 
the H2 emission in this galaxy is likely to be produced 
in PDRs. The small variations (a factor of 2) of the 
H2 S(l)/11.3 jam PAH ratio are probably due to differ- 
ences in the extinction or differences in the 11.3 /im PAH 
relative strength with respect to the total PAH luminos- 
ity from region to region. For the rest of the galaxies, the 
morphology of H2 S(l)/11.3 /im PAH ratio does not al- 
low us to determine whether the H2 is produced in PDRs 
or is excited by shock fronts or X-ray. 

8. VELOCITY FIELDS 

Figure [10] shows the velocity fields of the 
[Ne II] 12.81 /im and H2 S(l) emission lines in the 
SH module. They cover ~30"x30", which corresponds 
to physical scales of ~6-ll kpc, depending on the galaxy. 
Most of the [Ne II] 12.81 /im velocity fields are consistent 
with a rotating disk, although some of them show more 
complex morphologies (NGC 3110, NGC 3256, IC 694, 
NGC 6701 and NGC 7130). Only two of the galaxies 
(NGC 3256 and IC 694) with perturbed velocity fields 
are part of merger systems. NGC 3110 and NGC 7130 
are paired with other galaxies. However not all the 
galaxies in groups present disturbed velocity fields (e.g., 
NGC 7771, IC 4518E), at least on the physical scales 
(~kpc) probed by the IRS spectral maps. 

Figure [24] shows the velocity field of the interacting 
system Arp 299 (IC 694 and NGC 3690) rotated to the 
usual north up east to the left orientation for easy com- 
parison with other works. The [Ne II] 12.81 /im velocity 
field, which traces the ionized gas, agrees well with that 
measured from Ha (see Garcia-Marm et al. 2006). The 
peak-to-peak variation of the velocity field is '^250 km 
s~^ also comparable to that measured from Ha. The 
warm molecular gas (300 K < T < 1000 K) velocity 
field, traced by the H2 S(l) line, it is similar to the 
[Ne II]12.81 /im field (Figure [24]) , as both rotates in the 
same sense, although there are differences. The molec- 
ular hydrogen velocity field does not present the irregu- 
lar structure of the ionize d gas and the neutral gas (see 
iGarcia-Marm et"al]l2006f ). However, the H2 S(l) veloc- 
ity field resembles that produced by two disks rotating in 
opposite sense. This suggests that befor e the interactio n 
both galaxies were spirals (Augarde fc Lequeuxl Il985[ ). 
Moreover the H2 S(l) velocity field s hows evident simi- 
larities with the C O velocity field (see lAalto et al.iri997l : 
iCasoh et al1[T999[ ). 

Ha veloci t y fiel ds are available for NGC 3256 
(jLipari et al.l [2OOOI ). NGC 6701 and NGC 7771 
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Fig. 23. — Maps of the observed H2 S(l)/11.3 /im PAH ratio for selected galaxies. For reference we represent the 15 fim continuum 
contours. 
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Fig. 24. — Arp 299 (NGC 3690 and IC 694) velocity fields of the molecular gas traced by the H2 S(l) line (left) and the ionized gas traced 
by the [Ne II] 12.81 /im line (right). The velocities are referred to the the nucleus Bl of Arp 299. The black crosses mark the nuclei A and 
Bl of Arp 299. The scale represents 1 kpc. North is up and east to the left. 

to that of the [Ne II]12.81 /im. 

9. CONCLUSIONS 

We presented the analysis of Spitzer /IRS mapping ob- 
servations of 15 local LIRGs with extended Paa emission. 
We studied the spatial distribution of the mid-IR spectral 
features and compared the nuclear and integrated spec- 
tra. We calculated fine structure line and PAH ratios 
which trace the physical conditions in the star-forming 
regions. The main results are the following: 

• We used the 9.7 /im silicate feature strength vs. 
6.2 jam PAH equivalent width diagram to classify 
the activity of the LIRGs. There is a good agree- 
ment with the optical classification and this dia- 
gram. Most of the LIRGs populate the pure star- 



(lAlonso-Herrero et al.l 120093) and NGC 7591 
(jAmram et al.l 119941 ). The [Ne II] 12.81 /im velocity 
field of NGC 3256 appears perturbed due the merger 
nature of this galaxy. It is consistent with the Ha 
velocity field in shape and amplitude. The Ha velocity 
maps of NGC 7771 and NGC 6701 are slightly smaller 
than those presented here. NGC 7771 has and ordered 
velocity field compatible with a rotating disk whereas 
NGC 6701 presents a complex morphology both in 
Ha and mid-IR emission lines. The structure and the 
peak-to-peak velocity amplitude (~300 km s~^ for 
NGC 7771 and -150 km s"^ for NGC 6701) of the 
[Ne II] 12.81 /im velocity field are in good agreement 
with those of the Ha velocity field. As in the previous 
cases, the Ha velocity field of NGC 7591 is also similar 
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burst class and only the nuclei classified as Seyfert 
from optical spectroscopy appear in the diagram 
in the AGN/SB class. The integrated values show, 
in general, larger 6.2 fim PAH EW and shallower 
silicate absorption than the nuclear spectra. That 
is, the extended star formation partially masks the 
nuclear activity, resulting in a starburst-like inte- 
grated spectra. The silicate feature strength, nu- 
clear and integrated, of our sample of local LIRGs 
(S'si ^ —0.4 to — 0.9) is small compared to that 
found in ULIRGs (^si up to -4). 

We constructed maps of the spatial distribution 
of the [Ne III] 15.56 /im/[Ne II] 12.81 /im, which 
traces the radiation field hardness, and the 
[S III]18.71 /im/[NeIl]12.81 //m that varies with 
the electron density. In general, the minimum of 
both ratios is located at the nucleus. The nuclear 
[Ne III] 15.56 /xm/[Ne II] 12.81 fim ratio is below the 
expected range derived by star- formation models. 
A possible explanation is that due the high densi- 
ties in the nuclear regions this ratio is suppressed, 
including the possibilities that the most massive 
stars are either missing of buried in ultracompact 
Hll regions. Alternatively the star formation rate 
may have decreased rapidly in all nuclei over the 
last 10-20 Myr. 

We find a positive correlation between 
the [Ne III] 15.56 /im/[Ne II] 12.81 fim and 
[S III]18.71 /im/[Ne II]12.81 /im ratios for star 
forming regions in our sample of LIRGs 
and other starburst galaxies. On the other 
hand, AGNs, for a given [S III] 18.71 /im/ 
[Ne II] 12.81 /im ratio, show a systematically 
larger [Ne III] 15.56 /im/[Ne II] 12.81 /im ratio. 
The few starburst galaxies where we detect the 
[S IV] 10.51 /im emission line follow the correlation 
between the [S IV] 10.51 /im/[S III] 18.71 /im and 
[Ne III] 15.56 /im/[Ne II] 12.81 /im observed in Hll 
regions in other galaxies. We observe larger 
[S IV] 10.51 /im/[S III] 18.71 /im ratios for those 
galaxies harboring an AGN. 

We find that the 11.3 /im PAH emission is more ex- 
tended than that of the 5.5 /im continuum. How- 
ever, the ratio of the nuclear (2 kpc) 6.2 /im PAH 
emission with respect to the integrated emission 
is comparable to that of the 5.5 /im continuum 
in most cases. We find no correlation between 
the [Ne III] 15.56 /im/[Ne II] 12.81 /im ratio and the 
11.3 /im PAH EW, thus the effect of the radia- 
tion field hardness in the PAH emission, for the 
[Ne III] 15.56 /im/[Ne II] 12.81 /im ratio range in our 
sample of LIRGs, may not be important. We pro- 
pose that the 11.3 /im PAH/7.7 /im PAH ratio de- 
pends on the age of the stellar population. While 
the 7.7 /im PAH comes from young star- forming 
regions, the 11.3 /im PAH arises also from diffuse 
medium, thus the 11.3 /im PAH/7.7 /im PAH is 
lower in young Hll regions. We also found that 
large [Ne II] 12.81 /im/11.3 /im PAH ratios may in- 
dicate recent star formation. 

We explored the variation of the PAH ratios across 



the galaxies. In most cases these variations are 
real, that is, they are not due to extinction). In 
general the integrated 11.3 /im PAH/6.2 /im PAH 
ratios are larger than the nuclear values, proba- 
bly indicating that the integrated emission includes 
more neutral conditions. The nuclear and inte- 
grated 7.7 /im PAH/6.2 /im PAH ratios are almost 
constant, except in the nuclei of galaxies classi- 
fied as AGN where this ratio is higher. Since the 
6.2 /im PAH emission is more associated to small 
molecules, the increased 7.7 /im PAH/6.2 /im PAH 
ratio in AGN could be explained if the small PAH 
molecules are more easily destroyed in the harsh 
environments of active nuclei. 

• Using the H2 S(0) at 28.2 jam and the H2 S(l) 
at 17.0 /im lines integrated over 13"xl3" (3 to 5 
kpc) we estimated the mass of the warm (T^300 
K) molecular hydrogen. It ranges from 0.4 to 
2.6xlO^M0, and they are similar to those of 
ULIRGs, local starbursts and Seyfert galaxies. 
However these masses are lower limits since it is 
only included the H2 emission from the central few 
kpc of the galaxies and the total mass of warm 
molecular hydrogen is likely to be, at least, a fac- 
tor of 2 larger. 

• The similarity between the PAHs and molecular 
hydrogen morphologies suggests that the main ex- 
citation mechanism of the latter is the UV radiation 
too. However there are some regions with an excess 
of H2 emission with respect to the PAH emission, 
thus the other mechanisms should contribute no- 
ticeably to the H2 emission. Some of these regions 
are associated to interacting systems, where large 
scale shocks may also play a role in exciting the 
molecular hydrogen. 

• Despite the modest spectral resolution of the SH 
module, we show that useful velocity information 
can be obtained from the SH spectra. For most of 
the galaxies and on the physical scales probed by 
the IRS spectra (~kpc) the velocity fields are com- 
parable with that produced by a rotating disk. For 
the galaxies with available Ra velocity fields, the 
[Nell] 12. 81 /im velocity fields are in good agree- 
ment with those of the Ha, both in shape and peak- 
to-peak velocities. 
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TABLE 1 

Log of Spitzer/IRS Spectral Mapping Observations 







RA^ 


Dec^ 


Vhei^ Dist.'^ Scale 




Spect. 




Target 


IRAS Name 


(J2000.0) 


(J2000.0) 


(km s""*") (Mpc) (pc arcsec""*") 


T d 


class® 


AOR^ 



NGC 2369 


IRAS 


F07160-6215 


07 


16 


37. 


,75 


-62 20 37.5 


3196 


46. 


.0 


220 


11, 


.1 




17659392 


NGC 3110 


IRAS 


F10015-0614 


10 


04 


02. 


,12 


-06 28 29.1 


5014 


72. 


.5 


350 


11, 


.2 


Hii 


17659648 


NGC 3256 


IRAS 


F10257-4339 


10 


27 


51. 


,28 


-43 54 13.5 


2790 


40. 


,1 


190 


11, 


.6 


Hii 


17659904 


NGC 3690* 


IRAS 


F11257+5850 


11 


28 


30. 


,80 


+58 33 44.0 


3121 


44. 


.9 


220 


11, 


.8 


Sy2 


17660160 


IC 694* 






11 


28 


33. 


,80 


+58 33 46.5 














Hii 


17660416 


ESO320-G030 


IRAS 


F11506-3851 


11 


53 


11. 


,72 


-39 07 48.7 


3038 


43, 


.7 


210 


11, 


.2 


Hii 


17660672 


NGC5135 


IRAS 


F13229-2934 


13 


25 


44. 


,06 


-29 50 01.2 


4074 


58, 


.8 


280 


11, 


.2 


Sy2 


17660928 


IC4518E 


IRAS 


F14544-4255 


14 


57 


44. 


,80 


-43 07 53.0 


4584 


64, 


.5 


320 


11, 


.2 




21927936 


Zw049.057 


IRAS 


F15107+0724 


15 


13 


13. 


,12 


+07 13 31.7 


3858 


55. 


,6 


270 


11, 


.2 


Hii 


17661184 


IRAS17138-1017 


IRAS 


F17138-1017 


17 


16 


35. 


,80 


-10 20 39.4 


5161 


74. 


.7 


360 


11, 


.4 


Hii 


17661440 


IC4687 


IRAS 


F18093-5744 


18 


13 


39. 


,83 


-57 43 31.2 


5105 


73. 


.8 


360 


11, 


.5 


Hii 


17665024 


NGC6701 


IRAS 


F18425+6036 


18 


43 


12. 


,46 


+60 39 12.0 


3896 


56. 


.2 


270 


11, 


.0 


Composite 


21928448 


NGC7130 


IRAS 


F21453-3511 


21 


48 


19. 


,40 


-34 57 04.9 


4837 


69. 


.9 


340 


11, 


.4 


LINER/Sy2 


17661696 


IC5179 


IRAS 


F22132-3705 


22 


16 


09. 


,09 


-36 50 37.5 


3363 


48. 


.4 


230 


11, 


.2 


Hii 


17661952 


NGC7591 


IRAS 


F23157+0618 


23 


18 


16. 


,28 


+06 35 08.9 


4907 


70. 


.9 


340 


11, 


.0 


Composite 


21928960 


NGC7771 


IRAS 


F23488+1949 


23 


51 


24. 


,88 


+20 06 42.6 


4302 


62. 


.1 


300 


11, 


.4 


Hii 


21929472 



NGC 3690 and IC 694 are the two components, western and eastern respectively, of the interacting system Arp 299 
^ Coordinates of the nucleus of the galaxy. 

^ Heliocentric velocity from the IRS nuclear spectrum of each galaxy. 

Distance calculated using the measured heliocentric velocity and ou r cosm ology. 
^ Logarithm of t he L (8- 1000 fJ^m), as defined in Sanders & Mirabel (I1996D , in solar luminosities (3.826 X 10^^ W). It is calculated using the 
IRAS fluxes from [Sanders et al.| (I2003D and from Surace et al. (2004), the measured redshift, and our cosmology. 

® Classification of the nuclear activity from optica l spectroscopy. Composite indicates inte rmediate LINER/H ii class ification. The references for 
the nucle ar classification are given in lAlonso-H^rero et al.i (|2006j), except for NGC 3690 i Carcia-Marm et aI.ll2QQ6l ). NGC 6701 and NGC 7591 
(Alonso-H errero et al.ll2QQ9al) . 
^ Astronomical Observation Request key. 



TABLE 2 

Low RESOLUTION (SL) SPECTRA 



Nuclear 



Integrated 



PAH EW^ 



Name 



Size"" Ssi 6.2 /im 11.3 /im 



7.7 /iiii PAH 
6.2 /im PAH 



PAH EW^ 



11.3 /im PAH 
6.2 /im PAH 



Size'^ 



Ssi 



6.2 /im 



11.3 /im 



7.7 /im PAH 
6.2 /im PAH 



11.3 /im PAH 
6.2 /im PAH 



NGC2369 


1, 


.2 


-0.92 


0.43 


± 


0, 


.04 


0.33 


± 


0. 


.02 


1. 


.8 


± 


0. 


.0 


0.38 


± 


0. 


,01 


8.8x6.2 


-0.84 


0. 


.43 


± 


0. 


.08 


0.40 


± 


0. 


,03 


1, 


.8 


± 


0. 


,0 


0.55 


± 


0.01 


NGC3110 


2 


.0 


-0.48 


0.58 


± 


0. 


.11 


0.48 


± 


0. 


.05 


1. 


.7 




0. 


.0 


0.46 


± 


0, 


,01 


15.4x9.8 


-0.40 


0. 


.73 


± 


1. 


,02 


0.53 


± 


0. 


.05 


1, 


.6 


± 


0, 


,3 


0.49 


± 


0.06 


NGC3256 


1 


.1 


-0.48 


0.53 


± 


0. 


.02 


0.33 


± 


0. 


.01 


1. 


.7 


± 


0. 


.0 


0.46 


± 


0, 


,01 


7.6x3.8 


-0.52 


0. 


.59 


± 


0, 


,02 


0.43 


± 


0. 


.01 


1, 


.7 


± 


0, 


,0 


0.53 


± 


0.01 


NGC3690* 


1 


.2 


-0.80 


0.11 


± 


0. 


.01 


0.045 


± 


0. 


.001 


1. 


.4 


± 


0. 


.1 


0.21 


± 


0, 


,01 


8.3x5.7 


-0.96 


0. 


.29 


± 


0, 


,01 


0.21 


± 


0. 


.01 


1, 


.7 


± 


0, 


,1 


0.42 


± 


0.03 


IC694* 


1 


.2 


-1.9 


0.34 


± 


0. 


.01 


0.23 


± 


0. 


.01 


2. 


.1 


± 


0. 


.1 


0.28 


± 


0. 


,02 


8.3x5.7 


-0.96 


0. 


.29 


± 


0. 


,01 


0.21 


± 


0. 


.01 


1, 


.7 


± 


0. 


,1 


0.42 


± 


0.03 


ESO320-G030 


1 


.2 


-0.73 


0.49 


± 


0. 


.06 


0.50 


± 


0. 


.03 


1. 


.8 


± 


0. 


.0 


0.55 


± 


0. 


,01 


5.9x4.6 


-0.56 


0. 


.53 


± 


0. 


,10 


0.53 


± 


0. 


.04 


1, 


.8 


± 


0. 


,0 


0.63 


± 


0.01 


NGC5135 


1 


.6 


-0.47 


0.30 


± 


0. 


.02 


0.31 


± 


0. 


,02 


1. 


.7 


± 


0. 


.0 


0.58 


d= 


0. 


,01 


10x5.9 


-0.46 


0. 


.33 


± 


0. 


,05 


0.35 


± 


0. 


.03 


1, 


.7 


± 


0. 


,0 


0.73 


± 


0.02 


IC4518E 


1 


.8 


-0.82 


0.39 


± 


0. 


.28 


0.44 


± 


0. 


.17 


1. 


.8 


± 


0. 


.2 


0.47 


± 


0. 


,05 


6.4x6.4 


-0.65 


0. 


.62 


± 


0. 


,39 


0.53 


± 


0. 


.16 


1, 


.5 


± 


0. 


,1 


0.48 


± 


0.03 


ZW049.057 


1 


.5 


-1.2 


0.39 


± 


0. 


.15 


0.48 


± 


0. 


.16 


2. 


.1 


± 


0. 


.1 


0.49 


± 


0. 


,03 


4.9x3.8 


-1.2 


0. 


.42 


± 


0. 


,30 


0.47 


± 


0. 


.19 


2, 


.1 


± 


0, 


,2 


0.57 


± 


0.06 


IRAS17138-1017 


2 


.0 


-0.79 


0.56 


± 


0. 


.06 


0.40 


± 


0. 


.03 


1. 


.7 


± 


0. 


.0 


0.38 


± 


0. 


,01 


10x4.3 


-0.77 


0. 


.58 


± 


0. 


,10 


0.43 


± 


0. 


.04 


1, 


.7 


± 


0. 


,0 


0.48 


± 


0.01 


IC4687 


2 


.0 


-0.40 


0.66 


± 


0. 


.08 


0.55 


± 


0. 


.03 


1. 


.6 


± 


0. 


.0 


0.49 


± 


0. 


,01 


13x5.8 


-0.38 


0. 


.64 


± 


0. 


,11 


0.54 


± 


0. 


.04 


1, 


.7 


± 


0. 


,0 


0.57 


± 


0.01 


NGC6701 


1 


.5 


-0.54 


0.46 


± 


0. 


.07 


0.45 


± 


0, 


,05 


1. 


.7 


± 


0. 


.0 


0.58 


± 


0. 


,01 


13x6.5 


-0.40 


0. 


.53 


± 


0. 


,12 


0.52 


± 


0. 


.05 


1, 


.6 


± 


0, 


,0 


0.70 


± 


0.02 


NGC7130 


1 


.9 


-0.44 


0.25 


± 


0. 


.04 


0.23 


± 


0, 


,04 


1. 


.8 


± 


0. 


.1 


0.66 


± 


0. 


,02 


12x8.2 


-0.38 


0. 


.47 


± 


0. 


,10 


0.37 


± 


0. 


.05 


1, 


.5 


± 


0, 


,0 


0.65 


± 


0.02 


IC5179 


1 


.3 


-0.50 


0.57 


± 


0. 


.08 


0.46 


± 


0, 


,03 


1. 


.7 


± 


0. 


.0 


0.48 


± 


0. 


,01 


10x5.5 


-0.32 


0. 


.62 


± 


0, 


,08 


0.51 


± 


0. 


.03 


1, 


.6 


± 


0, 


,0 


0.54 


± 


0.01 


NGC7591 


1 


.9 


-0.85 


0.39 


± 


0. 


.12 


0.35 


± 


0, 


,06 


1. 


.9 


± 


0. 


.1 


0.52 


± 


0. 


,02 


16x8.8 


-0.47 


0. 


.63 


± 


0, 


,31 


0.43 


± 


0. 


.08 


1, 


.4 


± 


0, 


,1 


0.53 


± 


0.03 


NGC7771 


1 


.7 


-0.66 


0.46 


± 


0. 


.06 


0.38 


± 


0, 


,04 


1. 


.7 


± 


0. 


.0 


0.43 


± 


0, 


,01 


11x7.8 


-0.70 


0. 


.49 




0, 


,07 


0.41 


± 


0. 


.03 


1, 


.7 


± 


0, 


,0 


0.51 


± 


0.01 



NGC 3690 and IC 694 are part of the interacting system Arp 299. 
^ Projected square size of the aperture in kpc. The aperture size is 5.5" x5.5" for all the galaxies. 
^ Equivalent width units are /im. 

Approximated projected square size of the aperture in kpc. See Section l3.3l for details. 
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TABLE 3 

High resolution (SH) spectra. Integrated spectra 



Name 


Size^ 


[b ivj 
10.51/2m 


[Ne ii] 
12.81/im 


[Ne III] 
15.55/xm 


[S in] 
18.71/2m 


12.28/im 


H2 S(l) 
17.04/im 


11.3 fim PAH 


NGC2369 


9.7x 


:6.2 


<5.5 


194. 


± 11. 


15. 


± 3. 


44. 


± 4. 


19. ± 6. 


51. 


± 2. 


798. 


± 34. 


NGC3110 


15x 


8.4 


<6.9 


164. 


± 14. 


30. 


± 3. 


59. 


± 5. 


<5.1 


33. 


± 3. 


877. 


± 29. 


NGC3256 


6. Ox 


:3.8 


15. ± 3. 


1229. 


± 48. 


189. 


± 5. 


395. 


± 8. 


65. ± 6. 


124, 


. ± 3. 


2998. 


. ± 30. 


ESO320-G030 


6.3x 


:4.2 


<4.8 


152. 


± 7. 


12. 


± 2. 


47. 


± 5. 


13. ± 4. 


21. 


± 2. 


677. 


± 29. 


NGC5135 


9. Ox 


:5.0 


19. ± 2. 


141. 


± 12. 


66. 


± 3. 


39. 


± 3. 


14. ± 4. 


32. 


± 2. 


593. 


± 24. 


IC4518E 


8. Ox 


:6.4 


<1.6 


11. 


± 1. 


0.90 


± 0.32 


3.7 


± 1.0 


<1.0 


2.1 


± 0.5 


34. 


± 3. 


ZW049.057 


4.8x 


:4.3 


<4.3 


23. 


± 3. 


< 


2.2 


< 


1.4 


<6.1 


7.9 


± 1.0 


89. 


± 24. 


IRAS17138-1017 


11.5x5.8 


<3.9 


159. 


± 9. 


25. 


± 1. 


54. 


± 3. 


8.3 ± 3.7 


24. 


± 2. 


505. 


± 22. 


IC4687 


llx 


7.2 


10. ± 2. 


198. 


± 10. 


57. 


± 2. 


103. 


± 3. 


8.9 ± 5.5 


22. 


± 1. 


750. 


± 23. 


NGC6701 


8.1X 


:5.4 


<14. 


117. 


± 9. 


16. 


± 3. 


32. 


± 7. 


<8.8 


26. 


± 3. 


600. 


± 14. 


NGC7130 


16x 


:10 


<9.6 


158. 


± 14. 


36. 


± 4. 


39. 


± 7. 


<7.0 


29. 


± 5. 


643. 


± 29. 


IC5179 


lOx 


6.9 


<9.2 


290. 


± 17. 


28. 


± 3. 


91. 


± 5. 


19. ± 5. 


35. 


± 6. 


1334. 


. ± 29. 


NGC7591 


8.8x 


:6.8 


<5.8 


54. 


± 5. 


4.4 


± 2.5 


10. 


± 3. 


<3.6 


10. 


± 4. 


266. 


± 13. 


NGC7771 


14x6 


<12. 


148. 


± 11. 


<9.2 


41. 


d= 9. 


20. ± 10. 


41. 


± 4. 


675. 


± 23. 



Note. — Fluxes a re exp ressed in units of 10~ ^^Wm~^. The fluxes of Arp 299 (NGC 3690/IC 694) are given in 
lAlonso-Herrero et all (I2009br) . 

^ Approximated projected square size of the aperture in kpc. See Section l3.3l for details. 



TABLE 4 

High resolution (SH) spectra. Nuclear spectra 





[S IV] 


[Ne ii] 


[Ne III] 


[S III] 


H2 S(2) 


H2 S(l) 




Name 


Size"" 10.51/im 


12.81/im 


15.55/im 


18.71/im 


12.28/im 


17.04/im 


11.3 /im PAH 



NGC2369 


1.2 


<0.44 


52. 


± 3. 


3.4 ± 0.2 


8.6 ± 0.2 


3.1 ± 


0.2 


5.5 


± 0.2 


138. ± 3. 


NGC3110 


1.9 


<0.27 


26. 


± 1. 


2.4 ± 0.1 


6.2 ± 0.2 


1.3 ± 


0.3 


2.8 


± 0.1 


99. ± 2. 


NGC3256 


1.1 


1.8 ± 0.3 


260. 


± 14. 


27. ± 1. 


60. ± 1. 


8.0 ± 


0.8 


14. 


± 1. 


489. ± 5. 


ESO320-G030 


1.2 


<0.39 


41. 


± 2. 


3.3 ± 0.1 


11. ± 0.1 


2.8 ± 


0.3 


4.0 


± 0.1 


172. ± 3. 


NGC5135 


1.6 


8.0 ± 0.2 


51. 


± 3. 


23. ± 1. 


12. ± 0.1 


3.2 ± 


0.3 


4.9 


± 0.2 


166. ± 3. 


IC4518E 


1.8 


<0.082 


1.1 


± 0.1 


0.13 ± 0.02 


0.24 ± 0.08 


0.13 ± 


0.05 


0.34 


± 0.05 


5.4 ± 0.5 


ZW049.057 


1.5 


<0.26 


10. 


± 1. 


0.63 ± 0.07 


1.8 ± 0.2 


1.4 ± 


0.3 


2.2 


± 0.1 


32. ± 3. 


IRAS17138-1017 


2.0 


1.0 ± 0.2 


64. 


± 3. 


8.5 ± 0.2 


19. ± 1. 


2.1 ± 


0.6 


4.0 


± 0.4 


135. ± 3. 


IC4687 


2.0 


2.3 ± 0.2 


55. 


± 2. 


12. ± 0.1 


23. ± 0.1 


2.0 ± 


0.3 


3.2 


± 0.2 


189. ± 3. 


NGC6701 


1.5 


<0.49 


27. 


± 2. 


2.3 ± 0.1 


5.1 ± 0.2 


1.9 ± 


0.2 


3.4 


± 0.2 


103. ± 2. 


NGC7130 


1.9 


3.6 ± 0.4 


38. 


± 3. 


12. ± 1. 


7.4 ± 0.6 


1.7 ± 


0.8 


2.9 


± 0.3 


93. ± 3. 


IC5179 


1.3 


<0.39 


36. 


± 2. 


2.7 ± 0.1 


8.5 ± 0.3 


1.4 ± 


0.1 


2.6 


± 0.1 


120. ± 3. 


NGC7591 


1.9 


<0.24 


16. 


± 1. 


1.3 ± 0.1 


2.4 ± 0.2 


1.1 ± 


0.2 


1.8 


± 0.2 


48. ± 2. 


NGC7771 


1.7 


<0.38 


37. 


± 2. 


3.0 ± 0.1 


6.7 ± 0.4 


3.0 ± 


0.5 


4.8 


± 0.2 


115. ± 4. 



Note. — Fluxes are exp ressed in units of 10~^^ Wm~^. The fluxes of Arp 299 (NGC 3690/IC 694) are given in 
lAlonso-Herrero et all (l2009bD . 

^ Projected square size of the aperture in kpc. The aperture size is 5.5" X 5. 5" for all the galaxies. 



TABLE 5 
Nuclear flux/ Integrated flux 



Central 2.0 kpc/ Integrated flux Central 2.0 kpc/ Integrated flux 

5.5 /im 6.2 /im PAH 11.3 /im PAH 15 /im [Ne ii] 12.81 /im [Ne iii] 15.56 /im H2 S(l) 

Name continuum continuum 



NGC2369 


0.41 


± 


0.02 


0.41 


± 


0.01 


0.32 


± 


0.01 


0.39 


± 


0.02 


0.36 




0.02 


0.36 


± 


0.02 


0.20 


± 


0.02 


NGC3110 


0.13 


± 


0.05 


0.12 


± 


0.02 


0.10 


± 


0.01 


0.13 


± 


0.01 


0.07 




0.03 


0.12 


± 


0.02 


0.08 


± 


0.03 


NGC3256 


0.38 


± 


0.01 


0.36 


± 


0.01 


0.34 


± 


0.01 


0.39 


± 


0.01 


0.31 


± 


0.01 


0.41 


± 


0.01 


0.20 


± 


0.03 


ESO320-G030 


0.57 


± 


0.02 


0.54 


± 


0.01 


0.49 


± 


0.01 


0.48 


± 


0.01 


0.49 


± 


0.02 


0.40 


± 


0.01 


0.38 


± 


0.02 


NGC5135 


0.45 


± 


0.01 


0.42 


± 


0.01 


0.35 


± 


0.01 


0.39 


± 


0.01 


0.36 


± 


0.01 


0.28 


± 


0.01 


0.19 


± 


0.03 


IC4518E 


0.21 


± 


0.10 


0.18 


± 


0.10 


0.15 


± 


0.03 


0.14 


± 


0.02 








0.15 


± 


0.01 


0.19 


± 


0.06 


ZW049.057 


0.53 


± 


0.08 








0.50 


± 


0.02 


0.52 




0.02 








0.47 


± 


0.02 








IRAS17138-1017 


0.31 


± 


0.03 


0.31 




0.02 


0.25 


± 


0.01 


0.31 


± 


0.01 


0.27 




0.01 


0.19 


± 


0.01 


0.15 


± 


0.05 


IC4687 


0.24 


± 


0.03 


0.25 


± 


0.01 


0.23 


± 


0.01 


0.19 


± 


0.01 


0.17 


± 


0.01 


0.15 


± 


0.01 


0.13 


± 


0.05 


NGC6701 


0.25 


± 


0.03 


0.23 


± 


0.02 


0.20 


± 


0.02 


0.22 


± 


0.01 


0.16 


± 


0.02 


0.24 


± 


0.01 


0.15 


± 


0.03 


NGC7130 


0.28 


± 


0.02 


0.13 


± 


0.03 


0.14 


± 


0.02 


0.14 


± 


0.01 


0.21 


± 


0.03 


0.10 


± 


0.01 


0.09 


± 


0.05 


IC5179 


0.23 


± 


0.03 


0.20 


± 


0.01 


0.18 


± 


0.01 


0.19 


± 


0.01 


0.16 


± 


0.02 


0.11 


± 


0.01 


0.13 


± 


0.03 


NGC7591 


0.22 


± 


0.05 


0.14 




0.05 


0.13 


± 


0.02 


0.22 




0.01 


0.23 


± 


0.03 


0.23 


± 


0.01 


0.17 


± 


0.04 


NGC7771 


0.24 


± 


0.02 


0.28 


± 


0.02 


0.23 


± 


0.01 


0.24 


± 


0.01 


0.20 


± 


0.03 


0.23 


± 


0.01 


0.12 


± 


0.04 



Pereira-Santaella et al. 



TABLE 6 

Molecular hydrogen emission line fluxes 







H2 


S(2) 


H2 S(l) 


H2 S(0) 


Name 


Size^ 


12.27 i^m 


17.03 fim 


28.22 fim 


NGC2369 


2.9 


10. 


± 1. 


19. ± 2. 


3.5 ± 0.5 


NGC3110 


4.6 


4. 


± 1. 


10. ± 0.9 


1.6 ± 0.2 


NGC3256 


2.7 


34. 


± 3. 


61. ± 5. 


9. ± 3. 


NGC3690 


2.9 


11. 


± 4. 


22. ± 1. 


3.8 ± 0.6 


IC694 


2.9 


20. 


± 3. 


29. ± 1. 


5.7 ± 0.9 


ESO320-G030 


2.9 


8.1 


± 1.2 


13. ± 1. 


1.7 ± 0.7 


NGC5135 


3.9 


9.3 


± 1.4 


17. ± 1. 


1.3 ± 0.5 


IC4518E 


4.4 


0.3 


± 0.1 


1.1 ± 0.2 


1.0 ± 0.1 


ZW049.057 


3.7 


3.3 


± 0.8 


5.8 ± 0.4 


< 0.6 


IRAS17138-1017 


4.9 


6.3 


± 1.3 


14. ± 1. 


1.9 ± 0.5 


IC4687 


4.9 


6.5 


± 1.1 


12. ± 1. 


1.8 ± 0.4 


NGC6701 


3.7 


4.6 


± 0.8 


10. ± 1. 


1.9 ± 0.9 


NGC7130 


4.6 


4.2 


± 1.7 


9.9 ± 0.7 


1.4 ± 0.3 


IC5179 


3.2 


6.3 


± 0.8 


11.2 ± 0.8 


1.4 ± 0.2 


NGC7591 


4.6 


3.0 


± 1.0 


5.5 ± 0.6 


1.0 ± 0.3 


NGC7771 


4.1 


8.7 


± 1.2 


16. ± 1. 


3.5 ± 1.8 



Note. — Fluxes are expressed in units of 

Projected square size of the aperture in kpc. The 
aperture size is 13.4" x 13.4" for all the galaxies. 



TABLE 7 

Derived warm molecular hydrogen temperatures, column densities and masses. 





T(S(0)-S(1)) 


No 


Mo 


T(S(1)-S(2)) 


Ni 


Ml 


Name 


(K) 


(xlO^o cm-2) 


(xlO^ Mq) 


(K) 


(xlO^'O cm-2) 


(xlO^ Mq) 


NGC2369 


168 


8.8 


1.3 


348 


5.1 


0.73 


NGC3110 


178 


3.7 


1.3 


305 


2.4 


0.86 


NGC3256 


181 


19. 


2.1 


369 


12. 


1.3 


NGC3690 


172 


8.9 


1.2 


351 


5.3 


0.72 


IC694 


163 


15. 


2.0 


415 


7.8 


1.1 


ESO320-G030 


189 


3.4 


0.43 


390 


2.2 


0.28 


NGC5135 


236 


1.9 


0.44 


361 


1.5 


0.36 


IC4518E 


108 


9.2 


2.7 


252 


3.0 


0.87 


ZW049.057 


>205 


<1.0 


<0.22 


373 


0.74 


0.16 


IRAS17138-1017 


185 


4.1 


1.5 


324 


2.7 


1.0 


IC4687 


178 


4.1 


1.5 


361 


2.5 


0.92 


NGC6701 


168 


4.7 


1.0 


323 


2.8 


0.61 


NGC7130 


183 


3.0 


0.97 


320 


2.0 


0.65 


IC5179 


194 


2.6 


0.40 


369 


1.7 


0.27 


NGC7591 


167 


2.5 


0.85 


362 


1.4 


0.48 


NGC7771 


157 


10.0 


2.6 


363 


5.2 


1.3 



Note. — No and Mo are the column density and mass derived from the H2 lines S(0) and S(l). 
Ni and Mi correspond to the column density and mass derived from the S(l) and S(2) lines. 
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Fig. 1. — Spitzer /IRS SL spectral maps of the 5.5 /im continuum, the 6.2 /im PAH, the 7.7 /im PAH and the 11.3 /im PAH. The white 
cross marks the coordinates of the nucleus as listed in Table [T] The image orientation is indicated on the maps for each galaxy. The scale 
represents 1 kpc. The maps are shown in a square root scale. 
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Fig. 1. — Continued. 



Local Luminous Infrared Galaxies. Spatially resolved observations with IRS 
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Fig. 1. — Continued. 



26 



Pereira-Santaella et al. 



NGC2369 NGC3110 NGC3256 




arcsec arcsec arcsec 



Arp299 ESO320-G030 NGG5135 




arcsec arcsec arcsec 



10451 8E ZW049.057 IRAS17138-1017 




arcsec arcsec arcsec 



Fig. 2. — Spitzer/IRS SL spectral maps of the 11.3 /im PAH/6.2 /im PAH ratio. The 5.5 /im continuum contours, in a logarithm scale, 
are displayed to guide the eye. The image orientation is indicated on the maps for each galaxy. The scale represents 1 kpc. The ratio maps 
are shown in a linear scale. 
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Fig. 3. — Spitzer/IRS SL spectral maps of the 11.3 /im PAH/7.7 /im PAH ratio. The 5.5 fim continuum contours, in a logarithm scale, 
are displayed to guide the eye. The image orientation is indicated on the maps for each galaxy. The scale represents 1 kpc. The ratio maps 
are shown in a linear scale. 
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Fig. 4. — Spitzer/IRS SH spectral maps of the 11.3 /im PAH feature, the 15 fim continuum, the fine structure line emission of 



Local Luminous Infrared Galaxies. Spatially resolved observations with IRS 

NGC3256 



10 



Pa-a 



-10- 





|- Pnntikiiiim IRiirVi 



10- 



0- 




-10- UN 



10 



-10 



-10 10 

arcsec 



-10 10 

arcsec 



NGC3690 




lolPa 



-10-N^ 






10- 



[SIV]'l6.5|im' 



-10 





-10 

arcsec 



Fig. 4. — Continued. 
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Local Luminous Infrared Galaxies. Spatially resolved observations with IRS 
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32 



10- 



Pa- 




-10- V^N 



I I I I I I 



Pereira-Santaella et al. 
ZW049.057 



1 1 .3)im PAH 







UN 



Continuum 15)iin 



'UN 



- [Nell) 12.'8|im 



UN 



- [Nelli] is.slam 



UN 



10 



[SIM] 18.7)im 



IB 



-10- UN 



-10 



10 



- H2^(2) \ 2.3iirl\ 



14 



UN 



-10 



10 



nH7^(iyT70)!irrr 



UN 



-10 



10 



Pa-a 



I ' ' ' ' I 




-10-. /^r 



I I I I I I 



^Q^[SIV]'10.5)im 



0- 



I I I I I I 



I ' ^ ' ' I 



^Q^[SIII]'18.7)im 



0- 




I I I I I I 

-10 



IRAS17138-1017 

- n.'sjim pah' ' ' ' ' ' ' 




I I I I I I 



' [Nellj l'2.8|im' 




H2's'(2') 12!3|irT\ 




10 



I I I I I I 

-10 



10 




Fig. 4. — Continued. 



Local Luminous Infrared Galaxies. Spatially resolved observations with IRS 



34 



Pereira-Santaella et al. 
IC4687 



10 



Pa-d 



U 



-10- /^^L 

N ^ 



^Q^[SIV]'10.5| 



w 



10- 



-10- N 

, 






-10 10 
arcsec 



11.3|i^ 



"N E 



-10 




H2 'si2) l'2.'3|irT] 



~N E 








1 ' ' 

1 


^3 




T, , ,- 



-10 



-0[ 



--10 



-10 10 

arcsec 



NGC6701 



-10 

arcsec 



10 



10 



-10 



10 



-10 



'S(2)'i2'.3'|im' ' ' ' ' ' ' ' ' ' ' ^ 

h I I I I I I I I I I I I I I I I I I I 



-20 -10 10 20 




-20 -10 10 
arcsec 



Fig. 4. — Continued. 



Local Luminous Infrared Galaxies. Spatially resolved observations with IRS 



NGC7130 




arcsec arcsec 



Fig. 4. — Continued. 



36 




10L[SIII]18| 



n 




10- 



-io: 




Pereira-Santaella et al. 
NGC7591 




NGC7771 



PAH 




10 10 

arcsec 



l61 





-10 10 
arcsec 



-10 10 
arcsec 



10 10 
arcsec 



Fig. 4. — Continued. 



Local Luminous Infrared Galaxies. Spatially resolved observations with IRS 



37 



NGC2369 




-10 
arcsec 





arcsec 



NGC3110 




NGC3256 






arcsec 





arcsec 



NGC3690 






arcsec 





arcsec 



IC694 






arcsec 





arcsec 



Fig. 5. Spitzer/IRS SH spectral maps of the [S III]18.71 /im/[Ne Il]12.81 /xm ratio, the [Ne III]15.56 /im/[Ne Il]12.81 /im ratio and the 
[Ne II] 12.81 /im/11.3 fim PAH. The 15.0 /im continuum contours are displayed to guide the eye. The image orientation is indicated on the 
maps for each galaxy. The scale represents 1 kpc. The ratio maps are shown in a linear scale. 
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Fig. 8. — Maps of the 9.7 /im silicate feature strength as defined in Section \3A\ The 5.5 /im continuum contours are displayed to guide 
the eye. The image orientation is indicated on the maps for each galaxy. The scale represents 1 kpc. The maps are shown in a linear scale. 



42 



Pereira-Santaella et al. 




-20 -10 10 



91 0- 




^ 



NGC2369 
85 



13 

< 

-59? 
-131 
-203 

NGC3110 

269 




■ 175 



80 3 



-14 



arcsec 




-109 



NGC3256 
109 



68 



28 3 



-12 



-53 



121 



40 



-39 3 



-120 




arcsec 



-201 



224 



146 



- 68 3 




■ -87 



90 



48 



7 3 



-34 



-76 



10 



91 



-10 



^[Ne 


' 1 ' 
II] m 

ri 




1 ' ' 


' ' 1 
1 


1 








J 








1 




V 














E 

-, , , 1 , 


, , , 1 , , 


, , 1 





-10 



NGG3690 
145 





arcsec 



10 



81 



17 3 



-46 



-110 




172 



121 



71 3 



20 



-30 



Fig. 10. — Maps of the observed velocity field of [Ne II] 12.81 /imand H2 S(l) at 17.0 /im. The white cross marks the coordinates of the 
nucleus as listed in Tabled The image orientation is indicated on the maps for each galaxy. The scale represents 1 kpc. 



Local Luminous Infrared Galaxies. Spatially resolved observations with IRS 

IC694 






arcsec 



V? — 




arcsec 



ESO320-G030 
97 



46 



-5 3 



-56 



-108 




arcsec 




NGC5135 
80 



48 



16 3 



-14 



-46 





IC4518E 
221 



124 



27 3 



-69 



-166 




Fig. 10. — Continued. 



44 



Pereira-Santaella et al. 




ZW049.057 
134 



86 



39 3 



V? — 




-55 



IRAS17138-1017 
84 



42 



-41 



-83 



arcsec 




167 




IC4687 



10 



o 

n 



-10 



[Nell] 



^1 




I r 



9 3 



-69 



-148 



NGC6701 
124 



84 



43 3 



90 



59 



27 3 




-36 



78 



25 



-26 3 



-79 



arcsec 




-131 



137 



78 



19 I 



-40 



-99 



-20 -10 10 



20 



-37 




-20 -10 



10 20 



Fig. 10. — Continued. 



Local Luminous Infrared Galaxies. Spatially resolved observations with IRS 

NGC7130 




-20 -10 10 20 -20 -10 10 20 

arcsec arcsec 



NGC7591 




-20 -10 10 20 -20 -10 10 20 

arcsec arcsec 

Fig. 11. — Continued. 



46 



Pereira-Santaella et al. 



APPENDIX 

A. THE POINT SPREAD FUNCTION OF IRS SPECTRAL MAPPING OBSERVATIONS 

When comparing maps of lines or features, at different wavelengths, the variation of the Point Spread Function (PSF) 
size can be significant. The Spitzer /IRS PSF is much less studied than that of the other Spitzer instruments, although 
it is possible to generate synthetic IRS PSF using STinyTim program (Krist 2002). This is useful for staring obser- 
vations, but the spectral mapping observation mode introduces extra issues, such as the cube assembly from staring 
observations, wavelength calibration uncertainties, that cannot be easily taken into account using the modeled PSF. 
Because of that, a better way to characterize the IRS PSF of the spectral mapping; observation niode i s to use the 
mapping observations of calibration stars available in the Spitzer archive. Previously, iKospal et al.l (|2008f ) studied the 
factor that should be applied to correct for fiux losses when the slit is not centered properly on a point source. We 
extend this analysis to characterize the variati ons of the PSF. 

To do this, we used the same data sets as Ko spal et al.l (|2008f ). of spectral mapping observations of calibration stars. 
These observations were taken with the same instrument configuration used for the galaxies studied in this paper, but 
with much smaller steps between pointings, i.e. the calibration observations have a better spatial information of the 
PSF. Details of the observations are listed in Table [51 

After obtaining BCD files processed by Spitzer pipeline (versi on 15.3 for SH and SL modules, and version 17.2 for LH 
and LL modules), we built the spectral cubes with CUBISM (S mith et al.ll2007al ). The sky emission was subtracted 
using dedicated background observations. To take advantage of the better spatial information of the calibration 
observations, we set the pixel size in the final cube to be the size of the step in the perpendicular directiorO (see Table 
[8|). We did this for all IRS modules, however we only present the results for the SH, SLl and SL2 modules, because 
the other modules (LH, LLl and LL2), are heavily affected by fringing patterns and their PSF is not well reproduced 
by a Gaussian. 

In order to study the PSF core variations we fitted each plane of the cube using a 2-D Gaussian. We allowed different 
values of the Full Width Half Maximum (FWHM) for each axis of the Gaussian, but the axis orientation is fixed and 
aligned with the parallel and perpendicular directions. By doing this we obtained the FWHM and centroid of the PSF 
at every wavelength. These values are represented in Figure [25l 

The undulating pattern of the FWHM and the centroid in the modules SLl and SL2 is originated by several causes: 
(1) The spectral axis is not aligned with the detector, thus to extract the spectra is necessary to resample the original 
image. This resampling process produces this kind of periodic patterns. (2) The algorithm used to build the data 
cubes (see Smith et al.l l2007al for details) does not try to correct for this effect. In any case, these variations are small: 
<10 % for the FWHM and <0.2 arcsec for the centroid position. 

The PSF core is slightly asymmetric for all modules. The FWHM in the perpendicular direction is larger than the 
FWHM in parallel direction, which has an almost constant value. A possible explanation to this is the way CUBISM 
builds the cubes. It assumes that there is no spatial information in the dispersion (or perpendicul ar) direction and th e 
spatial information in that direction comes from different telescope pointings (for more details see iSmith et al.ll2007aD . 
The averaged values of the FWHM of each module are listed in Table [H 

The variation of the centroid position in SL modules is small. The la deviation around the average centroid position 
are 0.07 and 0.02 arcsec for the parallel and perpendicular direction respectively, that is a twentieth of the SL pixel 
size. These variations can be attributed to small telescope pointing errors and signal-to-noise limitation. The position 
of the centroid in the SH module shows variations that cannot be explained just by instrumental pointing error or a 
limited signal-to-noise. The l-cr deviation in the perpendicular direction is 0.03 arcsec, but in the parallel direction the 
variation is one order of magnitude larger (0.2 arcsec) and shows a clear wavelength dependency. The most probable 
cause for this, is an error in the position and/or angle of the extraction aperture. However this effect (0.6 arcsec) is 
small when compared with the SH pixel size (2.26 arcsec). 

We studied the variations of the diffraction rings in the SH module. We fitted the central row of the cube for each 
wavelength using a Gaussian for the PSF core and two Gaussians for the only diffraction ring detected. These two 
Gaussians have the same FWHM and they are symmetrically placed around the PSF core position. Figure [26l shows a 
fit example. The results of the fits are summarized in the Figure [27l The FWHM of the diffraction ring, that does not 
depend on the wavelength, is 3.6 ± 0.3 arcsec. Instead the position of the ring depends linearly on wavelength. The 
relative value of the diffraction ring peak respect to the PSF core peak, depends slightly on wavelength. At 10 /im it 
is 11% of the core peak value and at 19 /im it is 9% . 

When calculating feature ratios, the point spread function variations (FWHM and centroid), could originate artifacts. 
For the SL module these variations are less than one tenth of the pixel size and should not affect feature ratios. For the 
SH module the major source of uncertainty is the centroid position, but the difference is only one fourth of the pixel 
size. Since these variations are very small compared with the pixel size, they do not affect the calculated feature ratios. 
Diffraction rings can produce artifacts in the feature ratios too. They cause inaccurate ratio values 3 or 4 pixels away 
from the center of the point source in the SH module. For extended sources the situation is more complex because it 
is not possible to know where the diffraction ring will appear. However, comparing the spatial profiles of the features 
and the underlying continuum, as we do in this paper, helps determine if the ratio values are real or artifacts. 

In this Section we use the terms: "perpendicular" and "par- . . -, • r ^.i r^. j.- i 

n i„ m^ r ^ ^i t ^- i xi ^ i t pcrsiou axis and the spatial axis or the slit, respectively, 

allel direction, i hey refer to the direction along the spectral dis- 
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TABLE 8 

Log of IRS mapping calibration observations 



Step size perpendicular Step size parallel 



Module 


AOR 


Target 


(arcsec) 


(arcsec) 


SLl, SL2 


16295168 


HR 7341 


0.6 


15.5 


LLl, LL2 


16463104 


HR 6606 


1.7 


46.0 


SH 


16294912 


HR 6688 


0.8 


0.8 


SH 


16340224 


sky 






LH 


16101888 


HR 2491 


1.7 


1.7 


LH 


116088320 


sky 







TABLE 9 
IRS PSF Full Width Half maximum 





FWHM Perpendicular 


FWHM Parallel 


Module 


(arcsec) 


(arcsec) 


SH 


5.16 ± 0.08 


4.31 ± 0.37 


SLl 


3.76 ± 0.07 


3.27 ± 0.21 


SL2 


3.88 ± 0.04 


3.02 ± 0.17 



FWHM Parallel 



FWHM Perpendicular 



Center Parallel 



Center Perpendicular 




12 14 16 
Wavelength (|am) 



12 14 16 
Wavelength (|am) 



12 14 16 
Wavelength (|am) 



12 14 16 
Wavelength (|am) 




6.0 6.5 7.0 
Wavelength (|im 



6.0 6.5 7.0 
Wavelength (|am) 



6.0 6.5 7.0 
Wavelength (|am) 



Fig. 25. Full Width at Half Maximum (FWHM) and centroid position of the PSF of the modules SH, SLl and SL2. Parallel and 
perpendicular directions are represented. We have calculated these values by fitting a 2-D Gaussian. 
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SH 15.73|im 

i 




Fig. 26. — Upper panel: PSF of the SH module at 15.73 /im . Lower Panel: Fit of the PSF profile at the central row. The blue line is 
the Gaussian used to fit the PSF core. The green line is the fit of the diffraction ring. The red line is the sum of both components. 




Fig. 27. — Position of the diffraction ring in the SH module (black solid line). The position of the half-maximum of the diffraction ring 
is shown in gray. For comparison, the blue line is the half width at half maximum of the PSF core. The dashed black lines represent the 
pixel size of the SH module. 
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B. DESCRIPTION OF INDIVIDUAL OBJECTS 
B.l. NGC 2369 

The MIR emission of NGC 2369 is dominated by the resolved nucleus. The molecular hydrogen (II2 S(l) and H2 S(2)) 
emission appears more extended along the galaxy disk than the 11.3 jiia PAH emission. This causes a relatively large 
H2 S(l)/11.3 /im PAH ratio in the regions around the nucleus. The [Ne II] 12.81 /im/11.3 /im PAH ratio map reveals 
an Hll region ~17 arcsec (~3.4 kpc projected) away from the nucleus to the south. The nuclear silicate strength, 
5'si ~0.9, is slightly larger than the average silicate strength of our LIRGs sample, probably because this galaxy is 
seen almost edge-on. 

B.2. NGC 3110 

It is a face-on galaxy with numerous Hll regions along the spiral arms (|Hattori et ani2QQ4l ). The [Ne III] 15.56 /im 
emission of this galaxy is especially bright and it resembles that of the Ha emission. The [Ne III] 15.56 /im/ 
[Nell] 12. 81 /im ratio of the extranuclear Hll regions is ~5 times larger than the nuclear ratio. The deepest sili- 
cate absorption is located at the nucleus. 

B.3. NGC 3256 

It is galaxy merger with two nuclei. The southern nucleus, located 5 arcsec away from the north nucleus, is highly 
obscured, although it is visible in the 5.5 /im continuum map. It has one of the deepest silicate absorption feature 
in our sample, ^si ~ —1.4. The silicate feature of the rest of the galaxy is similar to that found in the rest of 
the galaxies of our sample. The PAH, continuum and atomic emission (except [S IV] 10.51 /im) peak at the north 
nucleus. The maximum of the [S IV] 10.51 /im emission is located ~5 arcsec to the east of the north nucleus, and it is 
probably related to a star-forining reg ion. This is the only extranuclear region (d>~l kpc), apart from region C in 
Arp 299 (|Alonso-Herrero et al.|[2Q"09b[ ). where we detect the [S IV] 10.51 /im line. The H2 emission arises mainly from 
the southern nucleus and the surrounding regions. 

B.4. ESQ 320-G030 

Due to the low spatial resolution of the IRS maps, we cannot distinguish any morphological detail in the emission 
maps this galaxy, although it is clearly spatially resolved. The maximum of the silicate feature absorption is located 
at the nucleus. 

B.5. NGC 5135 

This galaxy hosts a Sy2 nucleus surrounded by intense star formation (Bedregal et al.l [2009 1 and references therein), 
that we cannot separate due the limited spatial resolution of the IRS spectral maps. However the line ratios (e.g 
[Ne III] 15.56 /im/[Ne II] 12.81 /im ratio) show some asymmetries that may be associated with the truncated ring-like 
structure of this galaxy. At our scale (~0.5 kpc), we do not see any noticeable variations in the PAH emission due 
to the AGN activity. The silicate absorption feature is shallow (-0.47) and the maximum value is not located at the 
nucleus. 

B.6. IC 45I8E 

This galaxy is part of an interacting pair (IC 4518E/IC 4518W). Most of the mid infrared emission (PAH, fine 
structure line and 15 /im continuum) emission is dominated by an Hll region located ^7 arcsec to the south-east 
from the nucleus. The nuclear bright is similar to that of the southwestern H ll region emission only for the molecular 
hydrogen (H2 S(l)) and the 5.5 /im continuum. The deepest silicate absorption is found at the nucleus. 

B.7. Zw 049.057 

Due to the low signal-to-noise ratio of the observations of Zw 049.057, the maps of this galaxy contain very few 
pixels. The nuclear silicate absorption is, Ssi= -1.2, which is large when compared to that of the other LIRGs. 

B.8. IRAS 17138-1017 

The limited spatial resolution of the emission maps does not allow us to identify any morphological details in this 
galaxy. The maximum of the silica te absorption may not be located at the nucleus. Note however that the near- 
infrared (^ Alonso-Herrero et al.ll2006f ) and the mid-infrared (|Diaz-Santos et al.ll2008( ) morphologies do not allow us to 
determine the position of the nucleus. 

B.9. IC 4687 

IC 4687 is part of an interacting system (IC 4686/IC 4687). Only IC 4687 was observed in the spectral mapping 
mode. In the emission line map it is not possible to isolate any region of interest, however the emission is spatially 
resolved. We detect the [S IV] 10.51 /im line, that seems to be associated to the nuclear Hll regions seen in the Paa 
image. Two of these regions have slightly larger [Ne II] 12.81 /im/11.3 /im PAH ratio than the rest of the galaxy. The 
nuclear [Ne III]15.56 /im/[Ne II]12.81 /im ratio of this galaxy is somewhat larger than those of the other Hll galaxies 
in our sample. The maximum of the silicate feature is not located at the nucleus. 
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B.IO. NGC 6701 

This LIRG is classified as composite(H ll/LINER) from optical spectroscopy (|Alonso- Herrero et aIll2009aD . We can 
identify two extranuclear Hll regions symmetrically located ~10 arcsec away from the nucleus, one at each spiral arm. 
As is the case with the other extranuclear Hll regions, the molecular hydrogen of these two regions is relatively low 
when compared with the H2 nuclear emission. The silicate absorption is deeper in the nucleus than in the rest of the 
galaxy. 

B.ll. NGC 7130 

This is a face-on galaxy whose nuclear activity is classified as Sy/LINER. We detect 2 extranuclear Hll regions 
located in the spiral arms. The brightest is located to the north, ~10 arcsec from the nucleus. The other one is 
^14 arcsec to the south-east from the nucleus. The continuum and [Ne III] 15.56 /im emissions are concentrated at 
the nucleus. The PAH features and fine structure lines are bright in the spiral arms too. The minimum of the 
[Ne III] 15.56 /im/[Ne II] 12.81 /im ratio is not located at the nucleus, as for the majority of the LIRGs in our sample. 
We detect nuclear [S IV] 10.51 /im emission. Also, this is one of the few galaxies in the sample whose 11.3 /im PAH/ 
6.2 /im PAH ratio minimum is not located at the nucleus. These facts may be related to the active nature of the 
nucleus of this galaxy. 

B.12. IC 5179 

This is a face-on spiral gala xy classified as Hll-like from optical spectroscopy, and it shows numerous Hll regions 
distributed in the spiral arms (jAlonso- Herrero et al.|[2QQ6h . These Hll regions are not resolved in the Spitzer maps. 
They appear as a relatively uniform extended emission (PAH features and atomic lines) over ~8 kpc. The continuum 
emission is slightly more concentrated towards the nucleus. The patterns in the ratio maps are possibly related to the 
intricate distribution of the H ll regions. The nucleus is the most obscured region of this galaxy. 

B.13. NGC 7591 

The nuclear activity of NGC 7591 is classified as composite ( Alonso- Herrero et al.l '2009a'). Only in the low spectral 
resolution maps we can clearly distinguish the spiral arms of this face-on galaxy. The continuum emission is concen- 
trated at the nucleus, while we can see bright PAH and fine structure lines emission in the galaxy arms. The nucleus 
is the region with the deepest silicate absorption feature. 

B.14. NGC 7771 

The nucleus of this LIRG is classified as H ll-li ke from optical spectro scopy. It has a nuclear ring of star fo rmation 
clearly seen in the Paa images (Figure [H see also I Alonso- Herrero et al.l l2006 and I Alonso- Herrero et al.ir2009al ). but it 
is not resolved by the IRS observations. This galaxy is seen almost edge-on and we can clearly observe the spiral arms 
in the line emission maps. We detect an H ll region in the eastern arm, ~20 arcsec (6 kpc) away from the nucleus. The 
continuum is dominated by the nucleus. The H2 emission extends ~5 kpc in the west galaxy arm with apparently no 
bright fine structure/PAH emission associated with it. The silicate feature spatial distribution is somewhat irregular 
and we do not see any morphological detail. 
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